A 2 mm spectral line survey of the starburst galaxy NGC 253 
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ABSTRACT 

We present the first unbiased molecular line survey towards an extragalac- 
tic source, namely the nuclear region of the starburst galaxy NGC 253. The 
scan covers the frequency band from 129.1 to 175.2 GHz, i.e. most of the 2 mm 
atmospheric window. We identify 111 spectral features as transitions from 25 dif- 
ferent molecular species. Eight of which (three tentatively) are detected for the 
first time in the extragalactic interstellar medium. Among these newly detected 
species, we detected the rare isotopomers ^"^SO and HC^^O+. Tentative detec- 
tions of two deuterated species, DNC and N2D"'", are reported for the first time 
from a target beyond the Magellanic Clouds. Additionally, three hydrogen re- 
combination lines are identified, while no organic molecules larger than methanol 
are detected. Column densities and rotation temperatures are calculated for all 
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the species, including an upper limit to the ethanol abundance. A comparison 
of the chemical composition of the nuclear environment of NGC 253 with those 
of selected nearby galaxies demonstrates the chemical resemblance of IC 342 and 
NGC 4945 to that of NGC 253. On the other hand, the chemistries characteriz- 
ing NGC 253 and M82 are clearly different. We also present a comparison of the 
chemical composition of NGC 253 with those observed in Galactic prototypical 
sources. The chemistry of NGC 253 shows a striking similarity with the chem- 
istry observed toward the Galactic center molecular clouds, which are thought to 
be dominated by low-velocity shocks. This resemblance strongly suggests that 
the heating in the nuclear environment of NGC 253 is dominated by the same 
mechanism as that in the central region of the Milky Way. 

Subject headings: surveys — galaxies: abundances — galaxies: individual 
(NGC 253) — galaxies: ISM — galaxies: nuclei — galaxies: starburst — ra- 
dio lines: galaxies 

1. Introduction 

Our knowledge of the chemical composition of the interstellar medium (ISM) in the 
nuclei of external galaxies has so far been restricted to only ~ 30 molecular species (see 
Sect. 4.1), which is still far from the 129 molecules detected in the interstellar and circum- 
stellar medium within our Galaxy (Lovas 2004). The initial detection of most of these 129 
molecules has been the result of unbiased frequency scans towards specific Galactic sources, 
such as massive star-forming regions (e.g. Sgr B2(M), Cummins et al. 1986; OriMC-1, Lee 
et al. 2001), cold molecular clouds (e.g. TMC-1, Kaifu et al. 2004), and evolved stars (e.g. 
IRC-l-10216, Cernicharo et al. 2000). Spectral line scans towards Galactic sources focused 
mainly on the hot cores associated with massive star-forming regions such as SgrB2 and 
Orion-KL due to their complex chemistry as well as the brightness of their lines. Extra- 
galactic molecular spectroscopy has been so far limited to selecting the strongest features 
seen in the Galactic center and disk sources and observing them towards well selected extra- 
galactic targets. The most ambitious census of the molecular content of an external galaxy 
is the compilation of data toward the starburst galaxy NGC 4945 by Wang et al. (2004). 
Their multi-line study combines observations of 19 molecular species previously known in 
other extragalactic sources. 

Up to now unbiased molecular line surveys of extragalactic sources do not exist in spite 
of providing a powerful tool for determining the physical parameters of the molecular ISM. 
Observations of a large number of molecular lines with similar angular resolution allow us 
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to describe the chemical complexity of the source which provides fundamental information 
on the physical processes heating the medium in heavily obscured regions. Additionally, 
frequency scans lead to serendipitous detections of new species. 

The nuclear starburst of the galaxy NGC 253, the target of this survey, is one of the 
brightest and most prolific (in terms of molecular line detections) extragalactic molecular 
line sources. The Sculptor galaxy NGC 253 is an almost edgc-on barred spiral (i = 72° — 78°, 
Puche et al. 1991; Pence 1981) classified as type SAB(s)c (de Vaucouleurs et al. 1991) or 
Sc(s) (Sandage & Tammann 1987). At a distance of ~ 3Mpc (e.g. Mouhcine et al. 2005), 
NGC 253 is one of the nearest archetypes of nuclear starburst galaxies (Ricke et al. 1980). 
Its nuclear region contains one of the brightest extragalactic IRAS sources with a lOO/xm 
flux of 1860 Jy (Soifer et al. 1989). In fact, most of the overall IR luminosity of this galaxy 
(LiR = 2.1 10^°Lo,5Fi?iR = 3.6Moyr-^ Strickland et al. 2004) stems from the regions of 
intense massive star formation within its central few hundred parsecs. Violent massive star 
formation is also revealed by the high supernova rate of 0.05 — 0.3 yr"-*^ in the nuclear region 
of this galaxy (Mattila & Meikle 2001; Ulvestad & Antonucci 1997). 

The high nuclear star formation activity is driven by huge amounts of molecular gas in 
the central few hundred parsecs of NGC 253 (1.3 — 2.6 lO^M©, Canzian et al. 1988; Mauers- 

berger et al. 1996), likely powered by the bar structure first detected at near-infrared wave- 
lengths (Scoville et al. 1985; Forbes & DcPoy 1992). Position- velocity diagrams from high- 
resolution observations of the nuclear region of NGC 253 show that the molccTilar material 
seems to follow orbital motions around the dynamical center which arc interpreted as Xi and 
X2 orbits within the context of a bar-like potential (OH, Turner 1985; CO, Canzian et al. 
1988; Das et al. 2001; Paglione et al. 2004; HCN, Paglione et al. 1995; CS, Peng et al. 1996; 
SiO, Hi=^CO+, Garcia-BuriUo et al. 2000). 

Here we present the first unbiased molecular line survey carried out towards a source 
outside the Milky Way. Technically, this survey was driven by the availability of two new 
wideband (1 GHz) receivers at 2 mm and new 1 GHz wide filter bank spectrometers, together 
with the large collecting area of the IRAM 30m telescope. Species such as NO, NS, SO2, 
H2S and H2GS, detected for the first time towards a starburst environment, were previously 
reported by Martin et al. (2003, 2005). Additional new molecules together with a full analysis 
of the data are presented. 



2. Observations and results 
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Fig. 1. — Complete 2 mm spectral frequency scan toward the nuclear region of NGC 253. 
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The first molecular frequency scan of an external galaxy was carried out at 2 mm with 
the 30 m IRAM telescope at Pico Veleta, Spain, between 2001 and 2004. It covers ~86% 
of the observable 2mm atmospheric window, from 129.1 to 175.2 GHz (Fig. 1). At these 
frequencies, the telescope beamwidth ranges from 19" to 14". Fig. 2 shows the size of the 
30 m beam on top of an interferometric map of SiO emission (Garcia-Burillo et al. 2000). A 
i^-band image of NGC 253 (Engelbracht et al. 1998) is displayed in grey-scale to illustrate 
the size of the observed region relative to the galaxy. 
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Fig. 2. — Grey-scale, K-hand image of NGC 253 (left panel, Engelbracht et al. 1998) and 
high-resolution SiO emission (right panel, Garcia-Burillo et al. 2000). Beam sizes of the 
IRAM 30 m telescope for the extreme frequencies observed in the survey arc shown as circles 
on top of the SiO map. The continuous line shows the 19" beam at 129 GHz and the dashed 
line the 14" beam at 175 GHz. A cross indicates the nominal position of the 2 mm scan in 
this work and a dot the dynamical center as derived by Garcia-Burillo et al. (2000). 
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The observations were pointed at the position q;j2ooo = 00'^47™33!3, 5j2ooo = — 25°17'23". 
This position is ~ 6" south-east of the dynamical center of NGC 253. A pointing accuracy of 
~ 3" was achieved by measuring cross-scans on nearby continuum sources every ~ 2 hours. 
The observations were carried out in a wobbler-switching mode with a symmetrical beam 
throw of 4' in azimuth and a switching frequency of 0.5 Hz. 

Each of the two available SIS receivers working at 2 mm with orthogonal polarizations 
provides a 1 GHz bandwidth. The receivers were tuned to adjacent frequencies, with an 
overlap of 100 MHz, to cover an instantaneous bandwidth of 1.9 GHz for each frequency 
setup. As spectrometers wc used two 256 x 4MHz filterbanks, providing a velocity resolution 
between 7 and 9kms^^ at the observed frequencies. 

The rejection of the image band (the upper side band) was typically ~10 dB so that only 
the very few strongest lines were also detected from the image band. The image sideband 
rejection was calculated for each frequency setup (lower panel in Fig. 3) by measuring the 
difference in power between a hot and a cold load, first with full rejection of the upper side 
band and then with full rejection of the lower side band by using a Martin- Puplett interfer- 
ometer. The image gain is then computed as the ratio between these two measurements. We 
assume this gain ratio to be constant throughout the whole 1 GHz frequency band covered 
by each backend. 

The spectra were cahbrated with a standard dual load system. The temperature scale 
of the spectra is in Tmb obtained as 

Tub = {Fes / Bee)T^, (1) 

where the forward efficiency (-Feff) is 0.93 at 2 mm. The beam efficiency was calculated for 
each frequency with the Ruze function, 

Seff = 1.2eexp [-{iTTRa/X)% (2) 

using e = 0.69, Ra — 0.07 and A in mm. 

At the latitude of Pico Veleta Observatory, NGG253 is at elevations > 20° only for 
about four hours. A complete observing session with about 120 minutes effective observing 
time was spent on each frequency setup, out of which half the time was spent on source. 
The integration time for each frequency setup is shown in the central panel of Fig. 3. 

The observing mode combined with the stability of the system provided high quality 
basehnes. Only linear baselines were subtracted from the spectra. The rms of the residuals 
after subtracting the fitted line profiles on each frequency setup is shown in the upper panel 
of Fig. 3. The values are mostly of order of 2-4 mK in ~ 8 km s"^ wide channels, increasing 
up to ~6-10mK at the highest measured frequencies. 
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In the scanned 46 GHz wide band we have detected a total of 111 hnes, the weakest of 
which have an intensity of ~ 3 mK. The spectral density of detected lines is 2.4 features/GHz. 
Since the Doppler linewidth in NGC253 is about 200kms~^, or 100 MHz, we are not limited 
by line blending, but by a lack of sensitivity. For targeted searches with longer integration 
time it may still be possible to detect and identify lines at levels of < 1 mK. 

Fig. 4 shows the full 2 mm spectral line scan of NGC 253 with the line identifications 
superposed on the spectra. In Table 1 we are listing the parameters derived from Gaussian 
fits to the detected lines. We identify 25 molecular species and 3 hydrogen recombination 
lines. 
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Fig. 3. — Summary of the observational parameters of individual frequency setups. Open 
and filled symbols refer to each of the 2 mm SIS receivers. Upper panel: rms of the spec- 
trum in 7-9kms~^ wide channels in Tmb after subtracting a linear baseline; Central panel: 
Total integration time which includes the time spent on the reference position; Lower panel: 
Measured rejection of the image (upper) sideband. 
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Fig. 4. — 2 mm spectral line survey composite towards the nuclear region of NGC253. The 
resolution has been smoothed to four channels which is equivalent to velocity resolutions of 
35 to 27kms~'^. Line identifications are plotted on top of the spectra. Species labeled with 
i = molecule correspond to lines observed from the image band. 
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Table 1: Results from Gaussian fits to the observed lines 



Molecule 
Transition 


V 

(MHz) 


note 


J Tyi-B,Av 
mK km s 


1 —1 
kms 


1 —1 
km s 


Tmb 
mK 


SiO 3 — 2 v=0 


130268.6 




1580 ( 50) 


183 


76 


18 








1170 ( 50) 


260 


91 


12 


C2S llio - IO9 


131551.9 




640 ( 90) 


230 


151 


4.0 


HNCO 60,6 - 5o,5 


131885.7 




2400 (150) 


182 


72 


32 








2000 (180) 


284 


80 


24 


CH2NH 2i,i - li,o 


133272.1 


hf 


420 ( 40) 


187 


83*^ 


4.8 








450 ( 40) 


283 


83'' 


5.1 


OCS 11 — 10 


133785.9 




1220 (130) 


212 


200 


5.8 


S02 82,6 - 81,7 


134004.8 




300 ( 70) 


242 


91 


3.1 


H36o 


135286.0 


b s 


2800 (...) 


250* 


282* 


9.3 


±l2L^O 4i_4 — 3l,3 


135297.8 


b 


1040 ( 90) 


221 


158 


6.2 


SO2 5l,5 — 4o,4 


135696.0 




740 (160) 


245* 


140* 


4.2 


■^*SO 43 — 32 


135775.3 


b 


300 ( 60) 


180* 


73* 


3.8 






b 


200 ( 70) 


274* 


88* 


2.1 


HC3N 15-14 


136464.4 




2410 (130) 


183 


77 


29 








2000 (300) 


271 


85 


22 


CH3CCH 8fc - 7k 


136728.0 




3410 (110) 


255 


158 


20 


SO 43 — 32 


138178.5 




1400 (200) 


180 


73 


18 








1 -\ r\f\ i ^\r\f\\ 

1100 (200) 


274 


88 


12 


'^■^CS 3 — 2 


138739.3 




550 ( 70) 


188 


70*^ 


7.4 








390 ( 50) 


274 


70*^ 


5.2 


NH2CN 7i,7 - 61,6 


139032.0 




580 (110) 


280* 


155* 


3.5 


±l2i-^0 4l^3 — 3l^2 


ioy4oo.4 


D 


1 Q/i A ^1 
io4U i_10UJ 




I/O 


i .A 


NH2CN 7o,7 - 60,6 


139842.1 




490 (160) 


280* 


155* 


3.0 


NH2CN 7o,7 - 60,6 


139954.4 


m 


520 (160) 


279 


140 


3.5 


C2S lOii - 9io 


140180.7 




470 (130) 


250 


85 


5.2 


S02 62,4 — 61,5 


140306.1 




510 ( 60) 


248 


115 


4.1 


H2CO 2l,2 - ll,l 


140839.5 




4250 (140) 


186 


104*' 


38 








6050 (150) 


288 


104^ 


55 


C2S llii - lOlo 


142501.7 




430 ( 70) 


192 


85* 


4.7 


CH3OH 3i,3 - 2i,2 A+ 


143865.7 




1000 (130) 


225 


168 


5.6 


C2S 12ii — llio 


144244.8 




290 (120) 


210 


112 


2.4 


C^^s 3 _ 2 


144617.1 




3000 (300) 


188 


105 


27 








1500 (300) 


284 


86 


16 


C-C3H2 3l,2 — 22,1 


145089.5 


b s 


1750 (...) 


234* 


156* 


11 


CH3OH 3-2 


145103.2 


b m 


9400 (200) 


185 


115 


77 








5760 (190) 


294 


96 


57 


HC3N 16-15 


145560.9 


b s 


2100 (...) 


183* 


77* 


26 






b s 


1700 (...) 


271* 


85* 


19 


H2CO 2o,2 — lo,l 


145602.9 


b 


3400 (100) 


177 


90'' 


35 






b 


3300 (100) 


283 


90'' 


34 
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Table 1: (Cont.) 



Molecule 




v 


note 


J IVl D ^ 


Vlsr 


^^1/2 




Transition 




(MHz) 




mKkms""'^ 




kms""' 


mK 


OCS 12 - 11 




145946.8 




580 (170) 


212* 


210* 


2.6 






14636S 3 




1020 flSO") 


196 


139 


6.9 


SOo 4o O — 4l r> 




146605.5 




800 (170) 


241 


148 


5.1 


CS 3 - 2 




146969.0 


b 


11900 (200) 


185 


100 


111 








b 


13700 (200) 


288 


117 


110 


H35a 




147046.8 


b s 


2800 (...) 


250* 


282* 


9.3 


CHqCN 8t — 7t 




147174.5 


b m 


2040 (180) 


258 


165 


12 


HOCO+ 7n 7 - 6n fi 




149675.8 




410 (110) 


281 


68 


5.6 


NO 2 _ 1 n+ 

2 2 




150176.5 


hf 


3400 (500) 


244 


158 


10 


f—C iQ Ho A Til 

^ J ^ J ,U J- J -1- 




150436.5 


b s 


600 (...) 


234^ 


156* 


3.6 


H ol^tJ .zi 1 — 1 1 n 




150498.3 


b 


5300 (100) 


180 


99^ 


50 








b 


5100 (100) 


286 


99'' 


48 


NO 2 _ 1 rr- 

2 2 




150546.5 


hf b s 


3400 (...) 


244* 


158* 


10 


c-C iH'p 4n 4 — 3i 




150820.6 


b 


1150 (...) 


234'' 


156'' 


7.0 


c-C'-iH-j 4i zt - 3n ■? 




150851.9 


b 


3380 (110) 


234*^ 


156'' 


20 


C-0'^H'> 5i A — 5n c; 




151343.8 


b s 


165 ( ) 


234^ 


156* 


1.0 


oCJo 2o n — 2i 1 

■-'^-'Z '-'ZjU ^ i. , J. 




151378.6 


b 


590 (170) 


225 


158 


3.5 


DNC 2-1 




152609.7 




490 (160) 


225 


61 


7.5 


CH3CCH 9t. - 




153817.2 


b 


4300 (200) 


271 


157 


26 


HNCO 7o,7 - 60,6 




153865.0 




2660 (150) 


183 


72 


34 








b s 


2240 (...) 


284=^ 


80* 


26 


N2D+ 2-1 




154217.0 




570 (170) 


219 


120 


4.5 


HNCO 7i,6 -6i,B 




154414.7 




540 (110) 


283^ 


80* 


6.4 


HC3N 17-16 




154657.3 




1640 (120) 
1400 (110) 


191 
272 


72'' 
72'' 


21 
18 


C-C3H2 32,2 — 2i_l 




155518.2 




820 (180) 


205 


118 


6.6 


CH30H 60,6-6-1,6 


E 


157048.6 




540 (120) 


275 


103'' 


5.0 


CH3OH 5o.5 -5_l,5 


E 


157178.9 


b 


960 (150) 


242 


103'' 


8.8 


CH3OH Jo,fe — 


E 


157270.7 


b m 


2110 (160) 


183 


103'' 


19 


J = 1..4;A; = J 








5700 (200) 


290 


103'' 


52 


SO 34 - 23 




158971.7 




700 (200) 


180* 


92 


7.2 








b 


1000 (200) 


274* 


84 


11 


NH2CN 80.8 - 7o.7 




159814.6 




750 (180) 


283*' 


155'' 


4.5 


NH2CN 80.8 - 7o.7 




159942.7 


m 


700 (200) 


283'' 


158'' 


4.1 


H34a 




160211.5 




2800 (200) 


250 


282 


9.3 


NH2CN 81,7 - 7l,6 




161000.3 


b 


1430 (170) 


280 


155* 


8.6 


NS 1 - 1 / 




161297.2 


hf 


1800 (500) 


199 


224 


3.1 




161697.2 


hf 


2400 (400) 


249 


278 


3.4 


HC3N 18-17 




163753.4 




1500 (400) 
700 (300) 


185* 
272* 


73 
98 


20 
6.6 
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Table 1: (Cont.) 



Molecule 


1/ 


note 




^LSR 


^"1/2 




Transition 


(MHz) 




mK km 


km s^-'- 


km s^^ 


mK 


CH3OH Ji_fc_i — Jo,fc E 


165050.1 


b m 


1700 (400) 


169 


72 


22 


J = 1..3;fc = J 






4700 (400) 


273 


100 


45 


SO2 52,4 — 5l,5 


165144.6 


b s 


640 (...) 


248'^ 


115'^ 


5.2 


CH3OH 4i,3 - 4o,4 E 


165190.5 


b 


900 (200) 


291 


90^ 


9.1 


SO2 7ij — 60,6 


165225.4 


b s 


900 (...) 


248=- 


115^ 


7.4 


CH3CN 9fc - 8fc 


165568.9 


b m 


1200 (300) 


258*^ 


166» 


7.0 


H2S li,o — lo,i 


168762.7 




1520 (180) 


180 


100 


14 








2250 ( 80) 


275 


88 


24 


H2CS 5i,5 -4i,4 


169113.5 




970 (200) 


270 


165=^ 


5.5 


CH3OH 32,3 - 2l,3 


170060.6 




2400 (400) 


248 


178 


12 


2-1 


170322.7 




500 (300) 


206 


46 


11 








500 (200) 


285 


49 


10 


CH3CCH lOfc - 9fc 


170905.6 




5100 (500) 


275 


102 


47 


HOCO+ 80,8 - 7o,7 


171055.9 


b 


1700 (300) 


295 


66 


25 


SO 44 - 33 


172181.4 




1300 (300) 


180=^ 


87 


14 








900 (300) 


274=^ 


74 


11 


H13CN 2fe - Ifc 


172677.9 


hf 


3200 (200) 


181 


74 


40 








2400 (200) 


279 


66 


34 


HC3N 19-18 


172849.3 




2800 (400) 


184 


74 


35 








1800 (400) 


269 


75 


22 


H13CO+ 2 - 1 


173506.7 




860 ( 90) 


170 


eo'' 


14 








1450 (100) 


269 


60'' 


23 


SiO 4-3v=0 


173688.3 




1070 (110) 


184^ 


76'' 


13 








1000 (120) 


270=- 


91^ 


10 


C-C3H 3i,2 — 2i,i 


174086.1 


hf 


900 ( 200) 


180» 


4lb 


6.5 








1600 ( 500) 


280=- 


40 


12 


HNi3C 2-1 


174179.4 




1050 (170) 


271 


60 


16 


C2H 2-1 


174663.2 


hf 


30000 (6000) 


162 


108 


82 








27000 (6000) 


282 


104 


88 



Note. — (see Appendix A for details): 
b - Blended line. 

s - Synthetic Gaussian to isolate a given feature. See text for details on the parameters used for deriving 

these profiles for each species. 

hf - Hyperfine structure. Frequency and intensity refers to the main component of the group, 
m - Multi-transition line. Frequency refers to the main component of the group. 



"Parameter fixed in the Gaussian fit 

''Parameters forced to have the same value in the Gaussian fit 



-17- 



Table 2: Physical parameters derived for all detected species. Source averaged column den- 
sities, rotation temperature and abundances relative to H2. 



Molecule T^t [X]/[H2] " 









(cm ^] 




(K) 


(xlO-9) 


SiO 


180 km 


1 


5.0 


1.0) 


10^^ 


7.4(0.7) 


0.07 




270 km s" 


1 


3.6 


0.9) 


IO12 


8.7(1.2) 


0.05 


C2S 






1.4 


1.2) 


10^3 


24(9) 


0.2 


NH2CN 






1.2 


0.5) 


IOI3 


67(13) 


0.2 


CH3CN 






2.0 


0.6) 


1013 


9.6(0.7) 


0.3 


C-C3H 


130 km s- 


1 


8.1 


1.6) 


IOI2 


12 


0.1 




280 km s" 


1 


1.5 


0.3) 


1013 


12 


0.2 


C-C3H2 






3.0 


6.0) 


1013 


9(8) 


0.4 


NS 






4.3 


0.6) 


1013 


7.2(1.0) 


0.6 


HC3N 


180 km s" 


1 


1.8 


0.5) 


1013 


33(4) 


0.3 




270 km s" 


1 


1.9 


0.6) 


1013 


24(3) 


0.3 


HOCO+ 






2.6 


1.8) 


1013 


12 


0.4 


H2CS 






4.4 


4.5) 


1013 


11(4) 


0.6 


SO2 






5.0 


2.0) 


1013 


15(2) 


0.7 


H2S 


180 km s" 


1 


2.3 


1.1) 


1013 


12 


0.3 




275 km s" 


1 


3.5 


1.4) 


1013 


12 


0.5 


CH2NH 


190kms- 


1 


2.5 


1.1) 


1013 


12 


0.4 




300 km s" 


1 


2.8 


0.9) 


1013 


12 


0.4 


HNC 


270 kms- 


1 


7.2 


1.2) 


1013 


12 


1 


SO 


180 kms" 


1 


4.5 


3.3) 


1013 


40(24) 


0.7 




270 kms" 


1 


3.1 


2.8) 


1013 


34(24) 


0.5 


HNCO 


180 km s" 


1 


5.7 


2.7) 


1013 


23(6) 


0.8 




280 km s" 


1 


4.8 


2.5) 


1013 


23 


0.7 


HCO+ 


170kms- 


1 


3.9 


0.8) 


1013 


12 


0.6 




270 km s" 


1 


6.4 


1.2) 


1013 


12 


1 


H2CO 


180 kms" 


1 


6.9 


0.9) 


1013 


27(2) 


1 




285 km s- 


1 


9.0 


1.7) 


1013 


34(4) 


1 


OCS 






2.5 


0.3) 


1014 


17(2) 


4 


HON " 


180 km s" 


1 


1.8 


0.3) 


1014 


12 


2 




280 km s" 


1 


1.4 


0.2) 


1014 


12 


2 


CS "''^ 


180 kms" 


1 


2.0 


0.2) 


lO" 


9.7(0.4) 


3 




280 km 


1 


1.4 


0.2) 


10" 


10.0(0.2) 


2 


CH3CCH 






4.3 


0.2) 


1014 


63(20) 


6 


CH3OH 






8.3 


0.3) 


1014 


11.6(0.2) 


12 


C2H 


leOkms" 


1 


6.8 


1.1) 


1014 


12 


10 




280 km s" 


1 


5.6 


0.9) 


1014 


12 


8 


NO 






4.0 


0.4) 


lOis 


6(2) 


60 



"Errors derived from the rotation diagrams without beam filling factor considerations (see Sect. 3.3). 
''Assumed N{ll2) = e.riO^^cm-^. See Table 6. 

■^Calculated from the observed ^^C isotope transition with ^^C/^'^C ~ 40 from Henkel et al. (1993). 
''Using additional transitions measured by Martm et al. (2003, 2005). 
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3. Data analysis 

3.1. Line identification and fitting 

In order to identify and fit the observed features we have used the rest frequencies from 
the molecular hne catalogs of Lovas (1992, 2004) and Pickett et al. (1998). The existing 
molecular line surveys towards the Galactic center sources SgrB2(0H) (Cummins et al. 
1986; Turner 1989, 1991) and SgrB2(N, M and NW) (Nummclin ct al. 1998, 2000) were also 
used as additional information to estimate the expected intensity of the fainter identified 
lines not included in the Lovas (2004) catalog. 

One of the main difficulties in identifying lines from single dish observations of ex- 
tragalactic sources stems from their large width, typically > lOOkms"^. In addition, the 
emission lines from the nuclear region of NGC 253 show two velocity components, roughly at 
180 and 280 km s"^, with linewidths of ~ 100 and 110 km s^^ respectively (from CS, Martin 
et al. 2005). These components arise from the two main molecular lobes separated by 10" 
and located on both sides of the nucleus. As a consequence, some of the observed transitions 
will appear partially or totally blended. Table 1 shows the results of the Gaussian fits to all 
the molecular transitions identified towards the NGG 253 nuclear region. In the third column 
of Table 1 we have included a note to indicate whether the transition is affected by blending. 
Appendix A describes the fitting procedures used for the different types of blending found 
in our spectral line survey. 



3.2. Column density determination 

Prom the observed main beam brightness temperature of the measured molecular lines 
one can estimate column densities and rotation temperatures (Trot) for each species (see 
Appendix B.l for a detailed discussion). All the necessary spectroscopic information required 
to derive these parameters (i.e. A^i, u, g^, and Z in Eq. B6) were extracted or derived from 
the JPL catalog (Pickett et al. 1998). 

In order to derive Trot and to extrapolate the column densities in the observed states to 
a total column density for a given molecule, more than one transition has to be measured. 
In the case that only one transition was detected in the 2 mm scan, an educated guess of the 
rotation temperature was made (see Appendix C). 

Table 2 shows the derived total column densities and rotation temperatures for the 
detected species in NGC 253. Source averaged abundances are also given, where a source 
averaged molecular hydrogen column density A^(H2) = 6.710^^cm~^ has been assumed (see 
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Sect. 4.4.1). 
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Fig. 5. — Rotation diagrams of detected species. When differentiated, velocity components 
are represented as filled squares and solid lines (fSOkms^^) and filled triangles and dashed 
lines (270kms~^). Unfilled markers represent blended observed transitions whose intensities 
have been estimated by inter- or extrapolation from the diagrams (see Appendix A for 
details). The NS, NO and OCS diagrams include additional transitions observed by Martin 
et al. (2003, 2005). 
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3.3. Uncertainty in column densities and rotation temperatures 

The main source of uncertainty, when deriving the column densities and rotation tem- 
peratures, arises from the uncertainty of the extent of the emitting regions. 

The source averaged brightness temperature (Tb) can be estimated from the measured 
main beam brightness temperature as 

Tb = %l Tub (3) 

where, in the approximation of a Gaussian source distribution of size Og observed with a 
Gaussian beam of size O^, the beam filhng factor rj^f = 91/ [91 + 9l) accounts for the dilution 
effect due to the coupling between the source and the telescope beam. 

In order to correct for beam dilution in the nuclear region of NGC 253, Mauersberger et 
al. (2003) smoothed the high resolution interferometric GS map from Peng et al. (1996) to 
convert the three observed transitions of GS by Mauersberger & Henkel (1989) to a common 
angular resolution of 32". Since we know the intensities of the GS lines with the original 
resolution {Tmb) and those corrected to a 32" beam (Tj^g) we can derive the beam filling 
factor. Prom the expression 'q^f'TuB = VbjT^^ we estimate the extent of the emitting region 
to be Os = 23", 21" and 19" for the CS J = 2 - 1, J = 3 - 2 and J = 5 - 4 transitions 
respectively. In our analysis we will then consider an equivalent source size of 20" to convert 
our main beam brightness temperature to source averaged brightness temperature. 

The other main source of uncertainty stems from the assumption that all the observed 
species arise from the same volume. The only way to confirm this would be high-resolution 
observations of different transitions of several species. The interferometric map of the J = 
2 — 1 line of GS by Peng et al. (1996) shows that most of the emission of this molecule is 
concentrated within a radius of ~ 20". The GS emission is seen to be clumped into four large 
main molecular cloud complexes symmetrically distributed with respect to the dynamical 
center and roughly aligned with the molecular bar. A recent high-resolution multiline study 
carried out towards IG 342 (Meier & Turner 2005) shows clear differences in the distribution 
of several important molecular tracers. Towards NGC 253, however, high-resolution maps 
of molecular transitions requiring different excitation conditions such as HCN, CO, SiO, 
H13C0+ and NH3 (Paghone et al. 1995, 2004; Garcia-Burillo et al. 2000; Ott et al. 2005) 
show very similar distributions, in good agreement with that of GS. Therefore it is plausible 
that for NGG253, as long as the bulk of the emission arises from the area confined by the 
starburst, the extent of the emission of different molecules is similar. 

We now address the question of how an error in the assumed source size will affect 
our determination of physical parameters. As derived from Eq. Bl (see Appendix B.l), the 
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calculated total column density is affected by the filling factor through Tb- The upper panel 
in Fig 6 shows the variation of the inverse of the filling factor (i.e. the conversion factor from 
Tmb to Tb from Eq. 3) as a function of source size, normalized to the value for a source of 
20". This value represents, in a first approximation, the factor by which we should multiply 
our derived column densities if the source extent would be different from that assumed. We 
have considered the extreme cases of the beam sizes at 129 and 175 GHz. The grey shaded 
area corresponds to source sizes between 10" and 30" where we expect the source size to be 
confined. In case of a source size larger than 20" column densities would be overestimated 
by only ~ 20%, while for sizes down to 10", the column density would be underestimated 
by less than a factor of 2. Only in case of an extremely clumped and compact emission 
where source sizes were well below 10", column densities would be dramatically affected. 
We have also considered the effect of source size on the derived rotation temperature from 
rotation diagrams. We will assume a case in which we derive Tj-ot = 15K (for a 20" source) 
by observing two transitions of a given species at 129 and 175 GHz respectively. In the 
lower panel in Fig. 6 we have plotted the relative change in the derived Trot for those cases 
where the transitions would have upper levels with an energy difference of 60, 30 and 15 K. 
It is obvious that the more affected T^^t is that derived from the data with the smaller 
dynamical range in energies. Prom Fig. 6 we see that the largest expected uncertainty in the 
determination of T^ot is smaller than 20%. 
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10 20 30 

Source size (") 

Fig. 6. — Upper panel: Variation of the filling factor normalized to that corresponding to a 

source extent of 20" as a function of source size. The continuous line corresponds to the beam 
at 129 GHz and the dashed line to the beam at 175 GHz. The shaded area highlights values 
with source sizes between 10 and 30". Lower panel: Effect of the assumed source size on the 
rotation temperature determined from rotation diagrams relative to the assumed Trot = 15K 
for 6s = 20" with two transitions measured at 129 and 175 GHz where the involved levels are 
separated by 60 (continuous line), 30 (dashed fine) and 15 K (dotted line). 
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4. Discussion 

4.1. Extragalactic molecular census 

Table 3 lists the detection of all the species known to date in the extragalactic ISM 
in chronological order. During the last three decades, 29 molecules and 11 rare isotopic 
substitutions have been detected outside the Galaxy. As a result of the frequency scan 
presented in this paper (see also Martin et al. 2003, 2005), 8 new molecules and 2 isotope 
(shown in boldface in Table 3) have been added to the census of known extragalactic species, 
which represents a ~ 30% increase in the total number of identified molecules. While H2S 
and H2CS were detected towards the nearby Large Magellanic Cloud (Hcikkila ct al. 1999), 
their identification in our survey (see Martin et al. 2005) represents the first detections 
towards a starburst galaxy. Note that CH2NH, HOCO+ and C3H (Table 3) are tentatively 
detected and need confirmation by observations of additional transitions. In Table 3 the 
known extragalactic species not detected in the 2 mm band are shown in italics. 
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Table 3: Known extragalactic molecules (see Sect. 4.1) 



Molecule 


Publication 
Year 


Telescope 


Ref. 


OH 


1971 


OVRO 


1 


H2CO 


1974 


64 m 


2 


CO 


1975 


11m 


3 




1975 


11m 


4 


H2O 


1977 


100 m 


5 


HCN 


1977 


11 m 


6 


H2 


1978 


2.3 m 


7 


NH3 


1979 


100 m 


8 


HCO+ 


1979 


7m 


9 


CH 


1980 


64 m 


10 


cs 


1985 


7m 


11 


C3H2 


1986 


43 m 


12 


CH+ 


1987 


1.4m 


13 


CH3OH 


1987 


30 m 


14 


CN 


1988 


30 m 


15 


C2H 


1988 


30 m 


15 


HNC 


1988 


30 m 


15 


HC3N 


1988 


30 m 


15,16 


HNCO 


1989 


30 m 


17,18 




1989 


30 m 


19 


ci»o 


1991 


llm/30m 


20 


c"o 


1991 


llm/30m 


20 


so 


1991 


15 m 


21 


N2H+ 


1991 


30 m 


22 


SiO 


1991 


30 m 


22 




1991 


30 m 


22 




1991 


30 m 


22 




1991 


30 m 


22 


CH3CCH 


1991 


30 m 


23 


CH3CN 


1991 


30 m 


23 


13 CS 


1993 


30 m 


24 


ocs 


1995 


30 m 


25 


HCO 


1995 


12 m 


26 


DCO+ 


1996 


15 m 


27 


DON 


1996 


15 m 


27 




1999 


15 m 


28 


H2S 


1999 


15 m 


29 


H2CS 


1999 


15m/30m 


29,30 


CO+ 


2000 


30 m 


31 
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Table 3: (Cont.) 



Molecule 


Publication 


Telescoyje 


Ref. 




\ eai' 






SO2 


2003 


30 m 


32 


NO 


2003 


30 m 


32 


NS 


2003 


30 m 


32 




2003 


30 m 


32 


HOC+ 


2004 


30 m 


33 


C2S 


2006 


30 m 


34 


CH2NH 


2006 


30 m 


34 


NH2CN 


2006 


30 m 


34 


HOCO+ 


2006 


30 m 


34 


C3H 


2006 


30 m 


34 




2006 


30 m 


34 



Note. — Table adapted and updated from Mauersberger & Henkel (1993) 
In roman: molecules observed in this survey; In boldface: molecular species not previously identified outside 
the Milky Way; In italics: species not measured in this survey. 

Telescopes.- (OVRO) Owens Valley interferometer; (64m) Parkes 64m; (11m) NRAO Kitt Peak 11m; 
(100 m) Effelsberg 100 m; (2.3 m) Steward observatory 2.3 m; (7 m) AT&T Bell Labs 7m; (43 m) NRAO 
Green Bank 43 m; (1.4 m) ESO 1.4 m Coude; (30 m) IRAM Pico Veleta 30 m; (15 m) ESO SEST 15 m; 

References. — (1) Weliachew (1971); (2) Gardner & Whiteoak (1974); (3) Rickard et al. (1975); (4) 
Solomon & de Zafra (1975); (5) Churchwell et al. (1977); (6) Rickard et al. (1977); (7) Thompson et al. 
(1978); (8) Martin & Ho (1979); (9) Stark & Wolff (1979); (10) Whiteoak et al. (1980); (11) Henkel & Bally 
(1985); (12) Seaquist & Bell (1986); (13) Magain & Gillct (1987); (14) Henkel et al. (1987); (15) Henkel et al. 
(1988); (16) Mauersberger et al. (1990); (17) Nguyen- Q-Rieu et al. (1989); (18) Nguyen- Q-Rieu et al. (1991); 
(19) Mauersberger & Henkel (1989); (20) Sage et al. (1991); (21) Johansson (1991); (22) Mauersberger & 
Henkel (1991); (23) Mauersberger et al. (1991); (24) Henkel et al. (1993); (25) Mauersberger et al. (1995); 
(26) Sage & Ziurys (1995); (27) Chin et al. (1996); (28) Chin et al. (1999); (29) Heikkila et al. (1999); (30) 
Martin et al. (2005); (31) Fuente et al. (2000); (32) Martin et al. (2003); (33) Usero et al. (2004); (34) This 
work. 
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4.2. Deuterated molecules in NGC 253 

The observation of deuterated molecules in external galaxies is interesting for several 
reasons. The interstellar D/H ratio depends on the degree of processing of the gas, since D 
was mainly produced soon after the big bang while nuclear processes during stellar evolution 
rather destroy than form deuterium. A low D/H ratio may therefore indicate that the 
observed gas component has been recycled in stars, and a higher ratio could be a hint for 
recent infall of poorly processed material. It has been observed that D is strongly enriched 
(by a factor of > 10^ compared to the cosmic D/H ratio) in several molecules observable at 
cm- and mm-wavelengths. From model calculations, such a high degree of fractionation is 
expected if the molecules are being formed in a relatively cool (Tkin < 30 K) environment (e.g. 
Millar 2002, and references therein). The fact that a high enrichment of deuterated molecules 
is also observed in hot molecular cores can be explained by evaporation of molecules from 
grain mantles (e.g. Pineau des Forets et al. 1989). While a value of 3(±0.4) 10~^ has been 
estimated for the primordial D/H ratio (Buries 2002, and references therein), a lower ratio 
of 1.4(±0.24) 10~^ is found in the Galactic neighborhood, and variations are shown on larger 
scales (Vidal-Madjar 2002). In molecular clouds close to the center of the Milky Way, Jacq 
et al (1999) find DCN and DCO"'" enhancements which are compatible with a D/H an order 
of magnitude smaller than toward the local interstellar medium. This indicates that the gas 
in the Galactic center region is highly processed. The enhancement of deuterated molecules 
in the Large Magellanic Cloud, on the other hand, is compatible with the local value of the 
D/H ratio (Chin et al. 1996). No deuterated molecules have been detected so far toward 
the central regions of other galaxies. Mauersberger et al. (1995) give an upper limit for the 
DCN/HCN abundance ratio of 410"=^ toward NGC 253. 

Our survey includes the J = 2-1 lines of DCO+, DCN, DNC and N2D+. All other 2mm 
hnes of deuterated molecules are expected to be much weaker. Table 4 shows the frequencies 
of these hnes, the rms noise and hmits or values of the integrated intensities. DCO+ and DCN 
are not detected. There is a tentative feature at the frequency of DNC and a more significant 
feature at the line frequency of N2D''". At the frequencies of DNC J = 2 — 1 no line from 
another molecule is expected to be strong enough to explain the observed emission. Close to 
the frequencies of N2D+ J = 2 — 1 is the transition frequency of OC'^'^S J = 13 — 12 (at 154.242 
GHz). However, from our detections of OCS and our estimates of the ^^S/'^'^S ratio (Martin 
et al. 2005), this emission is expected to be an order of magnitude weaker than what is 
observed. Another potentially blending line is the 7i_6 — 61^5, Vq = 1 transition of vibrationally 
excited HNCO (at 154.228 GHz). This line has been detected toward the Galactic hot core 
GlO.47-1-0.03 with an intensity of about 1/5 that of the corresponding HNCO ground state 
transition, compatible with a kinetic temperature of about 380 K (Wyrowski et al. 1999). The 
feature we are observing here is about as strong as the HNCO 7i_6 ~ 61^5 ground vibrational 
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state transition (Table 1). If the feature at 154.228 GHz comes from vibrationally excited 
HNCO this would either indicate a higher kinetic temperature or maser emission. Both 
would be remarkable in view of the large volume observed in our beam toward NGC 253. 
While all this suggests that we have detected, for the first time, one or two deuterated 
molecules beyond the Magellanic Clouds, we will nevertheless regard the measured features 
as tentative. 
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Table 4: Parameters of selected deuterated species 



Molecule 


Transition 


1/ 


rms ^ 


/TMBd^;'^ 


[DX]/[HX] 




J- J' 


(MHz) 


(mK) 


(Kkms-^) 




DCO+ 


2 - 1 


144077.3 


1.9 


<0.30 


< 4 10-'^ " 


DCN 


2 - 1 


144827.9 


1.1 


<0.18 


< 1 10-2 " 


DNC 


2 - 1 


152609.7 


3.1 


0.49 


< 1 10-2 " 


N2D+ 


2-1 


154217.0 


2.8 


0.57 


< 1.410-3 ^ 



Note. — See Sect. 4.2 for details. 

"rms computed for 20kms~^ wide channels. 

^3ct limit to the integrated intensity assuming a linewidth of 150kms~^. DNC and N2D+ detections are 
tentative. 

■^Using the ^^C bearing species measured in this survey and assuming ^^C/^^C ~ 40 (Henkel et al. 1993). 
''Using the J = 1 — N2H+ transition from Mauersberger & Henkel (1991). 
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Table 4 gives limits to the [DX]/[HX] column density ratios for HCO+, HCN and HNC 
based on our limits for the J = 2 — 1 transitions of the deuterated isotopomers, the integrated 
intensities of the ^^C bearing isotopomers measured in this survey and the ^^C/^^C ratio of 
40 obtained by Henkel et al. (1993). We have assumed that the excitation conditions of 
the deuterated and main species are similar. Intensities were corrected to account for beam 
dilution assuming a source size of 20" (see also Sect. 3.3). For N2D+, we have only the 
N2H+ J — 1 — line for a comparison (Mauersberger & Henkel 1991). The value for the 
N2D+/N2H+ column density ratio given in Table 4 assumes optically thin emission, making 
the same assumptions on the excitation as they have been made for DCN in the appendix 
of Mauersberger et al. (1995). 

We obtain hmits to the [DX]/[HX] values of 4 10"^ for HCO+, of 1 10"^ for HCN and 
11 10~^ for HNC. These limits are two orders of magnitude above the local interstellar D/H 
ratio and an order of magnitude larger than the [DCN] /[HCN] ratios estimated for SgrB2 
(Jacq et al 1999), but an order of magnitude lower than the values determined by Chin et 
al. (1996) in the Large Magellanic Cloud. The limits are consistent with the high kinetic 
temperatures determined toward the molecular interstellar medium of NGC253. We can 
exclude the scenario found in many Galactic hot cores, where recent grain mantle evaporation 
leads to high D/H ratios that were preserved from a time when the temperature of the gas 
was much lower than now. 



4.3. Alcohols in NGC 253 

In the nuclear region of the Milky Way widespread large abundances of the methanol 
(CH3OH) and ethanol (C2H5OH) alcohols have been reported (Requena- Torres et al. 2005; 
Martin-Pintado et al. 2001). This fact has motivated a careful search for the presence of 
C2H5OH in the 2 mm scan. Table 5 shows the frequencies and energies of its 14 brightest 
2 mm transitions assuming Trot = 12K, as that derived for methanol (see Table 2). We 
have computed the rms of the noise for 20kms~^ velocity resolution whenever the transition 
appeared not to be contaminated by other spectral features. Prom this rms, we estimate 
a 3a detection limit to the integrated intensity of the transitions. With the limit to the 
detection of the 152 GHz transition we derive an upper limit to the ethanol column density 
of A^(C2H50H) < 1.8 10^^ cm^^. This limit to the column density can be improved if we add 
up all the transitions which appear not to be blended. We estimate a total column density 
limit of Ar(C2H50H) < 9.6 lO^^cm-^. 

The methanol to ethanol abundance ratio derived for NGC 253 is >10. This lower limit 
is still a factor of 2 below the ratios derived for the GC clouds by Requena- Torres et al. 
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(2005). Thus, further high sensitivity observations are required to estabhsh if the methanol 
to ethanol abundance ratio is similar to that found in the GC clouds, as expected if alcohols 
are also ejected from icy grain mantles by shocks in this galaxy (Charnley et al. 1995). 
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Table 5: Limits to the detection of C2H5QH 



V Eu rms * J TmbcIw ^ 

(MHz) (cm-i) (mK) (Kkms'^) 



129665. 


,7 


16.6 


3.0 


< 


1 


.440 


131103. 


,0 


16.6 


1.8 


< 


0, 


.864 


131502. 


,7 


20.1 


C2S 








132935. 


,1 


6.9 


baseline 








133323. 


,4 


16.5 


CH2NH 








142285. 


,0 


25.8 


1.8 


< 


0, 


.864 


144734. 


,0 


9.3 


C34g 








147427. 


,4 


21.0 


2.9 


< 


1, 


.392 


152656. 


,8 


9.3 


3.4 


< 


1, 


.632 


159781. 


,8 


12.2 


NH2CN 








160699. 


,0 


31.4 


i=H2S 








161609. 


,0 


26.1 


NS 








173391. 


,3 


12.2 


3.6 


< 


1, 


.728 


174232. 


,9 


15.7 


C-C3H 









"rms computed for 20kms ^ wide channels. If the transition is blended, the species is indicated. 
^3cr limit to the integrated intensity assuming an emission width of 150kms~^. 
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4.4. The nuclear environment in NGC 253 

The detection of new molecules in the extragalactic ISM, far from being a mere census 
of known species, provides a complete description of the chemical complexity within the 
nuclear environment of galaxies and a tool to establish the origin of this rich chemistry and 
the heating of the molecular gas. 

In the following subsections we will compare the observed abundances in NGC 253 with 
those available for five external galaxies (Sect. 4.4.1) and five prototypical sources within the 
Milky Way (Sect. 4.4.2). The comparison focuses mainly on the abundance ratios between 
the different molecular species compared to those observed toward NGC 253. This will 
minimize the systematic error in the absolute fractional abundances due to the H2 column 
density determinations. The relative abundances between molecules will provide us with a 
reliable chemical differentiation giving us insights into the dominant chemical mechanism 
acting in each source. 



4.4- 1- Comparison with other galaxies 

We have selected the nearby galaxies NGC 4945, M82 and IC342 where most of the 
known molecular species have also been observed, and Maffei2 and NGC 6946 for which 
fewer species have been detected. 

Table 6 presents the integrated intensities of the ^^CO J — 2 — 1 transition that has 
been used to estimate the H2 source averaged column density. We used the conversion 
A^(H2) = 3.3 lO^Ocm-^ 7(i3C02-i) (Mauersberger et al. 2003). 
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Table 6: Molecular hydrogen source averaged column density 



Source 


Transition 




^beam 






J -J' 


Kkms-^ 


(") 




NGC253 


2 - 1 


52.4*^ 


34 


6.71022 


NGC4945 


2 - 1 


81.2*^ 


23 


6.4 1022 


M82 


2 - 1 


29.4^ 


34 


3.81022 


IC 342 


2 - 1 


17.2^ 


34 


2.2 1022 


Maffei2 


i^CO 2 - 1 


13.9'^ 


34 


1.81022 


NGC 6946 


i^CO 2 - 1 


22.2'= 


21 


1.5 1022 



Note. — A 20" source size was assumed for all targets. 
''lOm HHT by Maucrsberger et al. (2003). 
''15 m SEST by Wang et al. (2004). 
^30 m IRAM by Israel & Baas (2001). 
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Table 7 shows the derived fractional abundances relative to H2 towards the selected 
galaxies. Because of the frequency coverage of our survey, not all the species listed in Ta- 
ble 3 could be measured. Therefore, in Table 7, we have also included data of additional 
key molecular species from other studies, namely NH3, N2H"'" and CN. H2 and CO are not 
included as they are tracers of the total mass of the molecular material content. Species such 
as OH and H2O have also been detected towards NGC 253 and many other galaxies. How- 
ever, H2O exhibits maser emission that does not allow us to determine directly the column 
densities, while OH shows extremely complex behavior involving a mixture of absorption, 
thermal and maser emission (e.g. Turner 1985; Ho et al. 1987; Prayer et al. 1998; Bradford 
et al. 1999; Henkel et al. 2004; Goicoechea et al. 2005). Unfortunately, there is not enough 
data available on CH, CH+, HCO, C0+ and HOC+ for this comparison. 

In order to achieve the highest homogeneity in the derived abundances, we have collected 
all the line profile parameters from the available observations in the literature. Prom these 
values we have made our own estimates of the column densities for each molecule and each 
source with the assumptions used for NGC 253. 

We have assumed a source size Og — 20" for all these sources (see e.g. Sect. 3.3 for 
NGC 253 and Wang et al. 2004 for NGC 4945). This is a fairly simple assumption for some 
of the targets and may lead to considerable errors in the column density determination. 
As this factor is also included in the source averaged column density determination of the 
molecular hydrogen, the error in the resulting fractional abundance will, however, cancel out 
to first order. 

The top panel in Fig. 7 shows a graphical representation of the data in Table 7. Molecu- 
lar species are ordered based on increasing abundances in NGC 253. We have also considered 
the logarithmic difference of abundances in each selected galaxy with respect to the abun- 
dances derived in this work for NGC 253 as AX = log^Q (X/Xngc253)- This variable, AX, 
is plotted in the bottom panel of Fig. 7. The difference in the F-axis between two species 
provides a measure of their abundance ratio relative to that in NGC 253. 

M 82 shows relative abundances quite similar to those of NGC 253 in many but not all 
of the observed species. Evident chemical differences between both nuclear starbursts were 
already analyzed (Mauersberger & Henkel 1993; Takano et al. 1995, 2002; Wang et al. 2004), 
and have been interpreted as a difference in the evolutionary stage of both starbursts. Com- 
pared to that in NGC 253, the nuclear starburst in M 82 represents an evolved state where 
PDRs dominate the heating of the molecular material (Garcia-Burillo et al. 2002). Prom 
our comparison, the high relative abundance of C3H2, a common FDR tracer, is particu- 
larly outstanding. On the other hand we note remarkably low abundances of CH3OH and 
SiO, which are thought to be ejected from grain mantles into the gas-phase (Charnley et al. 
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1995; Martm-Pintado et al. 1992; Requena-Torres et al. 2005). This suggests that shocks 
do not play the dominant role within the nuclear region of M82. Molecules such as HNCO, 
apparently associated with shock chemistry (Zinchenko et al. 2000), also show low relative 
abundances compared to those in NGC 253. Furthermore, the nuclear environment of M 82 
is characterized by systematically lower relative abundances of CH3CN, NH3 and N2H"'". It 
is also interesting to note how other species used as PDR tracers, namely C^^S and C2H 
(Meier & Turner 2005), do not show high relative abundances in M82 compared to the other 
observed species. Moreover, they are a factor of ~ 2 — 3 less abundant than most other 
the species. If we compare the M82 relative abundances with those of the Orion Bar in 
Sect. 4.4.2, we find that all the species with low abundances in M82 show the same behavior 
in the Orion Bar. This comparison strongly supports the idea of M 82 as an evolved nuclear 
starburst largely dominated by PDRs. 

In NGC 4945, relative abundances of all the detected molecular species resemble within 
a factor of 3 those of NGC 253. The non-dctcction of SiO indicates a clear undcrabundance, 
similar to what is seen towards M82. On the other hand, unlike in M82, CH3OH shows 
an abundance similar to that in NGC 253. This resemblance between the abundances in 
NGC 4945 and NGC 253, altogether with excitation considerations and a detailed study of 
species such as CN, HCN, and HNC, has been interpreted as an intermediate evolutionary 
state, between NGC 253 and M82 (Wang et al. 2004). Prom our comparison, the nuclear 
starburst of NGC 4945 shows much closer resemblance to NGC 253 than to M82. 

IC 342 shows high relative abundances of HC3N and H^^CO^, while molecules such as 
C2H and N2H"'" have lower abundances as compared with NGC 253. Many of the available 
observations for IC 342 (marked by an asterisk in Table 7) were obtained near the D position 
as labeled by Downes et al. (1992), ~ 15" north of the main nuclear star formation region. 
The abundance of C2H given in Table 7 has been taken from Meier & Turner (2005) towards 
the D position. Its relative abundance is low compared to that in NGC 253. Even taking the 
C2H abundance at the position where this molecule peaks (i.e. a value an order of magnitude 
higher towards the central position, Meier & Turner 2005), the relative abundances of this 
species with respect to the rest of molecules would be similar to that in NGC 253. High 
resolution maps from Meier & Turner (2005) clearly differentiate a region dominated by UV 
radiation from the central nuclear cluster and the region where the chemistry is dominated 
by shocks produced by collisions of clouds in a barred potential. Given the position towards 
which many of the spectra were taken, away from the dynamical center, the molecular 
emission of IC342 should resemble more that of NGC 253 than that of M82. In fact IC342 
does not show the relatively low HNCO abundance seen in M 82 but does show a similar low 
relative abundance of C2H. Still, there is a large number of molecules not yet observed which 
are critical to fully evaluate the chemistry of IC342, such as SiO, SO, OCS and CH3OH. 
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As far as Maffei2 and NGC 694.6 are concerned, the available molecular information 
is not sufficient to establish a clear chemical differentiation from the other sources. What 
can be noted is a remarkably high relative abundance of HNCO in Maffei2 compared to the 
other galaxies and a low relative abundance of H2CO in NGC 6946. Both galaxies also show 
low relative abundances of N2H"'". 
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Table 7: Molecular fractional abundances toward NGC 253 and other galaxies 



Molecule 


NGC 


NGC 


M82 


IC 


Mairei2 


NGC 




253 


4945 




342 




6946 


HNi^C 


-10.6 


-9.5 


-9.5 








H13CO+ 


-10.4 


-10.0 


-9.9 


-9.5* 






H13CN 


-9.9 


-9.7 


< -9.9 


< -9.4* 






SiO 


-9.9 


< -9.9 


< -9.9 


< -9.3* 


< -9.1 




CH3CN 


-9.5 




< -9.7 










-9.4 


-9.4 


-9.3 


-9.5* 




< -9.0 


C-C3H2 


-9.3 


-8.9 


-8.1 








HC3N 


-9.2 


-8.8 


-8.7 


-8.6 






N2H+ 


-9.2 


-9.2 


-9.4 


-9.6* 


-9.7 


-9.6 


HNC 


-9.0 


-8.6 


-8.8 


-8.7 


-8.7 


-8.8 


SO 


-8.9 


-8.7 


< -8.5 


< -8.7 






HCO+ 


-8.8 


-8.4 


-8.4 


-8.8 


-8.6 


-8.5 


HNCO 


-8.8 


-8.4 


< -8.8 


-8.6 


-8.3 




H2CO 


-8.6 


-8.1 


-8.2 


-9.4 


-8.8 


-9.1 


OCS 


-8.4 


< -7.5 


-7.9 


< -8.6 






CN 


-8.3 


-7.7 


-8.2 


-8.3* 






HON 


-8.3 


-8.3 


-8.4 


-8.5* 


-8.3 


-8.5 


CS 


-8.2 


-8.4 


-8.2 


-8.4* 


-8.7 


-8.7 


CH3CCH 


-8.2 


-8.0 


-7.7 








CH3OH 


-7.9 


-7.4 


< -8.3 


< -7.8 


-7.7 


-7.9 


C2H 


-7.7 


-7.3 


-7.6 


< -8.4 






NH3 


-7.2 




-8.4 


-7.3 


-7.1 





*Observations tciken ~ 15" away from the central position (see Sect. 4.4.1 for details). 

References. — NGC 253, this paper, Martin et al. (2003)[NS,NO], Martin et al. (2005) [OCS], Mauersberger et al. 
(2003)[NH3], Hcnkcl ct al. (1988)[CN], Mauersberger & Henkel (1991)[N2H+]; 

NGC 4945, Wang ct al. (2004), Henkel et al. (1994) [SiO]; 

M82, Mauersberger & Henkel (1989)[CS,C-^'*S], Mauersberger & Henkel (1991)[N2H+,SiO,Hl-'^CO+], Mauersberger et 
al. (1991)[C3H2,CH3CN,CH3CCH], Mauersberger et al. (1995)[OCS], Hiittcmeister et al. (1997)[H2CO,CH30H], Henkel 
et al. (1988)[CN,HC3N,C2H], Henkel et al. (1998)[HCO+,HCN,Hi3CN], Petuchowski & Bennett (1992)[SO], Nguycn-Q- 
Rieu et al. (1991)[HNCO], Nguyen-Q-Rieu et al. (1992)[HCN,HCO+,Hi3cO+], Wild (1990)[HCN,HCO+], Weifiet al. 
(2001)[NH3],Hiittemeister et al. (1995)[HNC]; 

IC342, Mauersberger & Henkel (1989)[CS], Mauersberger & Henkel (1991)[N2H+,SiO,Hi3cO+], Mauersberger et al. 
(1995)[C'^*S,OCS,HC3N,H2CO], Henkel et al. (1988)[CN,HNC,HC3N], Huttemeister et al. (1995)[HNC], Hiittemeister et al. 
(1997)[H2CO,CH30H], Petuchowski & Bennett (1992)[SO], Wild (1990)[HCN,HCO+], Nguyen-Q-Ricu ct al. (1991)[HNCO], 
Nguyen-Q-Rieu et al. (1992)[HCN,HC0+,H1^C0+], Henkel et al. (1998)[HCN, H^^^CN], Meier & Turner (2005)[C2H], Mauers- 
berger et al. (2003) [NH3]; 

Mafrei2, Mauersberger & Henkel (1989)[CS], Mauersberger & Henkel (1991)[SiO,N2H+], Huttemeister et al. (1997)[H2CO, 
CH3OH], Nguyen-Q-Rieu et al. (1991)[HNCO], Nguyen-Q-Rieu et al. (1992) [HCO+, HON], Henkel et al. (2000); Mauersberger 

et al. (2003) [NH3]; 

NGC 6946, Huttemeister et al. (1997)[H2CO, CH3OH], Mauersberger & Henkel (1989)[CS, C^-^S], Nguyen-Q-Rieu et al. 
(1992)[HCO+, HON], Mauersberger & Henkel (1991)[N2H+]. 
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Fig. 7. — Molecular abundances (X) of selected extragalactic sources compared to those 
of NGC253. Arrows represent upper limits. A continuous line connects the abundances 
in NGC253. Top panel: Logarithmic fractional abundances relative to H2 as presented 
in Table 7. Bottom panel: Logarithmic abundances for each species normalized to that 
measured in NGC 253. 
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In order to quantify this comparison, Table 8 shows the mean and rms of AX for 
each source for two cases. The first only uses the species detected in both sources. The 
number of molecules included in the calculations is also indicated. The second also uses 
upper limits to the abundances derived for the comparison sources. To take these limits 
into account, we followed two criteria: upper limits higher than the abundances measured 
towards NGC 253 were rejected as they do not provide useful information for the comparison; 
for the calculations, the other observed limits, lower than the abundances in NGC 253, 
where divided by two as an estimated true value. The second assumption may overestimate 
some abundances, but any other constraint might result in biased guesses. While the mean 
value, AX, only gives an the overall absolute difference in the abundances of each source 
as compared with NGC 253, the rms of AX provides us with a fairly reliable measurement 
of differences in the abundance pattern that defines the chemistry of each object. M82, as 
pointed out before, clearly differs chemically from NGC 253, showing the highest rms values. 
NGC 4945 and IC342 are much closer to NGC 253. 
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Table 8: Statistical comparison of abundances of selected galaxies with those in NGC 253. 

Without limits With limits 



Source 


AX 


rms(AX) 


num. 


AX 


rms(AX) 


num. 


NGC 4945 


0.30 


0.29 


18 


0.27 


0.32 


19 


M82 


0.19 


0.51 


15 


0.05 


0.54 


19 


IC 342 


0.08 


0.34 


12 


0.00 


0.42 


13 


Maffei 2 


0.02 


0.31 


9 


0.02 


0.31 


9 


NGC 6946 


-0.20 


0.30 


7 


-0.20 


0.30 


7 



Note. — A.X = log^o (-''^/-''^NGC253)- See last paragraph in Sect. 4.4.1 for details. 
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4 ■4- 2- Comparison with prototype Galactic sources 

We have selected five Galactic sources, namely SgrB2(N), SgrB2(0H), TMC-1, L134N, 
and the Orion Bar, which are considered to be prototypes of their respective chemistry within 
the Galaxy. All these sources have been the target of multi-transition studies and some of 
them of frequency surveys. Table 9 shows the abundances of all the molecules observed 
in the NGC253 survey compared to those of these Galactic sources. In order to achieve 
consistency among ratios, we have tried, when available, to use values from line surveys or 
multi-line studies for each source. 

In the same way as in the extragalactic comparison, the top panel in Fig. 8 shows a 
graphical representation of the data in Table 9, where molecular species are ordered based 
on increasing NGC253 abundances. The bottom panel in Fig. 8 shows the logarithmic 
abundance of each molecule normalized to the abundance measured in NGC 253. 

SgrB2(N) is believed to represent the prototype of a hot molecular core chemistry as- 
sociated with massive star formation near the nucleus of the Milky Way. A large number of 
molecular species have been identified towards this source which is the brightest molecular 
line emitter within the Galaxy (Friedel et al. 2004, and references therein). There are clear 
differences between the abundances found towards SgrB2(N) and those in NGC 253. A set 
of molecules, namely CH3CN, HC3N, NS, H2CS, SO2, CH2NH, SO, CH3OH and NO, show 
abundances relative to the other species ~ 2 — 3 orders of magnitude larger than what we 
observe towards NGC 253. Almost all the other species observed towards SgrB2(N) have 
smaller abundances than in NGC 253, but their relative abundances are, within an order 
of magnitude, similar to those in NGC 253. Exceptions are HOCO+ and C2H that show 
lower relative abundances by factors ~30 and 100, respectively, in the SgrB2(N) hot core. 
Moreover, the rotation temperatures derived from the observed molecules in NGC 253 are 
low compared to the typical temperatures of > 70 K observed towards hot cores. Even if hot 
cores associated with massive protostars are present in the nuclear environment of NGC 253, 
they clearly do not dominate the molecular emission. While the angular size of the SgrB2 
molecular cloud complex is ~ 18' (Scoville et al. 1975), the emission of some large complex 
molecules towards SgrB2(N) appears concentrated within a region of ~ 5", corresponding 
to 0.2 pc at a distance of 8.5 kpc (Kerr, F.J. & Lynden-Bell, D. 1986). At the distance of 
NGC 253, such a source would have a diameter of < 0'.'02, and because of beam dilution, its 
detection is certainly out of reach for the 30 m telescope. A large number of such sources 
would be needed in order to be able to observe traces of their complex chemistry. 

The Orion Bar has been selected as the prototype of photodissociation regions (PDRs). 
A large fraction of the molecular gas in the Milky Way and in external galaxies is expected to 
be affected by PDRs (HoUenbach & Tielens 1997). In PDRs, molecular emission originates 
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from the surface layers of interstellar molecular clouds exposed to intense far-ultraviolet 
photons from nearby young OB stars. The non-detection of molecules such as CH3CN 
and HNCO towards the Orion PDR clearly indicates a low relative abundance of these 
species with respect to other molecules, which is in contrast with the observed NGC253 
abundances. Similarly, the relative abundances of SiO, SO2, CH3OH, C2H and NO appear 
to be ~ 1 — 2 orders of magnitude lower in the Orion bar than in NGC 253. High abundances 
of hydrocarbon radicals, such as C2H and C3H2, are observed towards PDRs (Fosse et al. 
2000; Fuente et al. 1993, 2003; Pety et al. 2005). In this context, it is surprising to see 
how C2H and C3H2, commonly used as PDR tracers, show relatively low abundances in 
the Orion PDR when compared not only with NGC 253 but also with the prototypical 
Galactic sources except the SgrB2(N) hot core. The observations of the Orion Bar have 
been taken from the position of maximum H2 column density (Jansen et al. 1995; Fuente 
et al. 1996), which does not match the position of highest C3H2 abundance (Fuente et al. 
2003). However, the contribution of the position of the C3H2 emission peak to the total 
column density is smaller than that of the H2 column density peak. Therefore, it is expected 
that the molecular abundance of the H2 peak dominates the emission of PDRs, also in 
external galaxies. Anyhow, even if we consider the position of the larger C3H2 abundance 
(up to one order of magnitude, Fuente et al. 2003), the relative abundance would still be close 
to that found in NGC 253. From this comparison as well as that to M82 (see Sect. 4.4.1) 
it is clear that the chemistry of the nuclear environment of NGC 253 is not dominated by 
photodissociation. 

Quiescent cold dark clouds are represented by TMC-1 and L134N. The molecular com- 
position at sites of low mass star formation is dominated by gas-phase ion-molecule chemistry 
due to a lack of embedded luminous sources. In our comparison, both dark cloud complexes 
have in common low relative abundances of SiO and, to a lesser degree, of CH3OH. The 
remaining species behave quite differently in each source. SO25 H2S, HNCO and NO, show 
low abundances in TMC-1 while C2S and C3H2 appear to have abundances larger than in 
NGC 253 by up to one order of magnitude. On the other hand, L134N shows a high SO 
abundance and relatively low abundances of HC3N and CH3CCH. The abundance pattern 
defined by these two dark clouds clearly does not follow the relative abundances found to- 
wards NGC 253. In addition, even the rotation temperatures derived for NGC 253, which in 
general are smaller than Tkin, are higher than the typical kinetic temperature (Tkm ~ 10 K) 
measured towards dark clouds. 

SgrB2(0H), located at the southern end of the SgrB2 molecular envelope, is free from 
the emission steming from the hot cores in this cloud. This position in the SgrB giant 
molecular cloud is taken as the prototype of Galactic center (GC) molecular cloud complexes 
(Martin-Pintado et al. 1997). The H2 column density of this source has been estimated to be 
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roughly A'"(H2) = 5 10^^ (Nummelin et al. 2000). The relative abundances of all the species 
measured towards this Galactic center molecular cloud complex appear to closely follow the 
abundance pattern in NGC 253 within a factor ~ 5. The only species clearly underabundant 
in SgrB2(0H) when compared to NGC 253 is C-C3H but this molecule is only tentatively 
detected towards NGC 253. Thus, the ratios between any of the observed species are similar 
to those in NGC 253 within one order of magnitude. 
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Table 9: Comparison of NGC253 fractional molecular abundances with those in selected 
Galactic sources 



Molecule 


NGC 


Sgr 


Sgr 


TMC-1 


L134N 


Orion 




253 


B2(N) 


B2(OH) 






Bar 




-10.6 


-11.0 










h1'^CO+ 


— 10.4 


— 11.4 








—10.3 


SiO 


-9.9 


-10.7 




< -11.6 


< —11.4 


—10.3 


NH2CN 


-9.7 


-10.1 


-10.0 








C2S 


-9.7 




-9.6 


-8.1 


-9.2 




CH3CN 


—9.5 


—6.7 


—9.4 


—9.0 


< —9.0 


< —10.3 


C-C3H 


-9.5 


-10.5 


< -10.9 


-9.3 








O A 

— y.4 


— lU.O 


n 7 

— y. ( 








C34S 


-9.4 


-10.2 








-9.0 


C-C3H2 


-9.3 


-10.5 


-9.8 


-8.0 


-8.7 


-9.7 


HC3N 


-9.2 


-7.5 


-9.0 


-8.2 


-9.7 




NS 


-9.2 


-7.0 




-9.1 


-9.5 




H2CS 


-9.2 


-6.8 


-8.7 


-8.5 


-9.2 




SO2 


-9.1 


-6.6 


-8.7 


< -9.0 


-8.4 


-9.9 


CH2NH 


-9.1 


-7.0 


-9.2 








H2S 


-9.1 


-9.9 




< -9.3 


-9.1 


-8.2 


HNC 


-9.0 






-7.7 


-8.2 


-9.0 


SO 


-8.9 


-6.9 


-8.7 


-8.3 


-7.7 


-8.0 


HCO+ 


-8.8 






-8.1 


-8.1 


-8.5 


HNCO 


-8.8 


-9.2 


-8.4 


-9.7 




< -10.8 


H2CO 


-8.6 


-9.3 


-8.6 


-7.7 


-7.7 


-8.2 


ocs 


-8.4 


-8.6 


-8.3 


-8.7 


-8.7 




HCN 


-8.3 






-7.7 


-8.4 


-8.3 


OS 


-8.2 






-8.0 


-9.0 


-7.6 


CH3CCH 


-8.2 


-8.4 


-8.8 


-8.2 


< -8.9 




CH3OH 


-7.9 


-5.8 


-7.3 


-8.7 


-8.5 


-9.0 


C2H 


-7.7 


-9.7 




-7.1 


< -7.3 


-8.7 


NO 


-7.2 


-6.0 




< -7.5 


-7.2 


-8.6 



References. — NGC253: this paper and Martin et al. (2003) [NS.NO], Martin et al. (2005)[OCS]; 
SgrB2N: Nummelin et al. (2000) and Minh et al. (1991)[H2S]. N{B.2) = S.OlO^^cm-^ from Nummelin et al. (2000); 
SgrB2(0H): Cummins et al. (1986) and Turner (1991). Ar(H2) = 5 10^^ cm-^ (see text for details on this estimate); 
TMC-1 and L134N: Ohishi et al. (1992), McGonagle et al. (1994)[NS], Ziurys et al. (1989)[SiO]. iV(H2) = 1.0 10^2 cm-^ from 
Ohishi et al. (1992); 

Orion Bar. Jansen et al. (1995), Puente et al. (2003)[C3H2], Schilke et al. (2001)[SiO]; N(H.2) = 6.5 10^2 cm-^ from Jansen et 
al. (1995). 
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Fig. 8. — Molecular abundances (X) of selected Galactic sources compared to those of 
NGC253. Arrows represent upper limits. A continuous line connects the abundances in 
NGC 253. Top panel: Logarithmic fractional abundances relative to H2 as shown in Table 9. 
Bottom panel: Logarithmic abundances for each species normalized to that measured in 
NGC 253. 
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In analogy to our statistical comparison of molecular abundances in prototypical nearby 
galaxies shown in Table 8, Table 10 presents a comparison between NGC 253 and the selected 
Galactic sources. To calculate the parameters, the available abundances in Table 9 have 
been taken into account following the criteria explained at the end of Sect. 4.4.1. Given the 
uncertainty of its detection, C-C3H was not included. The computed rms of the logarithmic 
differences shows the envelope of the SgrB2 complex as the one most closely resembling 
NGC 253. 
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Table 10: Statistical comparison of abundances of selected Galactic sotut:cs with those in 
NGC253. 

t 111 limits With limits 



Source 


AX 


rms(AX) 


num. 


AX 


rms(AX) 


num. 


SgrB2(N) 


0.40 


1.50 


22 


0.40 


1.50 


22 


SgrB2(0H) 


0.07 


0.40 


16 


-0.04 


0.58 


17 


TMC-1 


0.47 


0.66 


18 


0.26 


0.85 


21 


L134N 


0.18 


0.58 


16 


-0.00 


0.76 


18 


OrionBar 


-0.07 


0.71 


14 


-0.28 


0.85 


17 



Note. — 



See Table 8 and Sect. 4.4.2 for details. 
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The agreement in relative abundances between NGC 253 and a typical Galactic center 
cloud is a remarkable and clear indication that the heating and chemistry of the bulk of the 
ISM is dominated by the same processes. 

The molecular envelope around the star forming region Sgr B2 is formed by a relatively 
dense and warm material (riHa = 2 lO^cm"^ and Tkin = 40 — 60 K, de Vicente et al. 1997). The 
large abundances of NH3, SiO and C2H5OH observed in the gas phase of the Galactic center 
molecular clouds are claimed to be the result of grain chemistry and subsequent ejection 
onto the gas phase due to the disruption of grains by low velocity shocks (v < 20kms^^, 
Flower et al. 1995; Mart in- Pint ado ct al. 1997, 2001). In the Galactic center, the possible 
origin for these shocks is so far unclear. It has been claimed that they could be produced 
by the large-scale shocks associated with cloud-cloud collisions due to the orbital motion 
of molecular clouds in a barred potential (Hasegawa et al. 1994; Sato et al. 2000) and the 
interaction of expanding bubbles due to supernova events and/or strong stellar winds from 
Wolf-Rayet stars (Sofue 1990; de Vicente et al. 1997; Martin-Pintado et al. 1999) associated 
with a burst of star formation that occurred 10^ years ago (Rodriguez-Fernandez et al. 2004; 
Rodriguez- Fernandez & Martin-Pintado 2005). 



5. Conclusions 

1. We present the first unbiased molecular line survey of an cxtragalactic source. The 
survey covers the 2 mm atmospheric window from 129.1 to 175.2 GHz toward the inner 
200pc of NGC 253. A total of 111 features arc identified as transitions of 25 different 
molecular species, eight out of which (three tentatively) have been detected for the first 
time outside the Milky Way. The rare isotopes, ^^SO and HC^^O"*", were also detected 
for the first time in an cxtragalactic source. Additionally, three hydrogen recombination 
hues and, tentatively, two deuterated species are identified, N2D+ and DNC. The 
origin of the observed features of both N2D+ and DNC is still unclear and deserves 
further investigation. If real, these would represent the first deuterated molecules 
observed beyond the Magellanic Clouds. Column densities and rotation temperatures 
have been determined for each species under Local Thermodynamic Equilibrium (LTE) 
conditions. 

2. As a result of this survey and by adding existing data from three molecules, namely 
NH3, N2II+ and CN, we obtain the most complete description to date of the chem- 
ical complexity within the nuclear few hundred parsecs of a starburst galaxy. We 
have performed a comparison between the chemistry of NGC 253 and those of five 
other outstanding nearby galaxies. This comparison clearly shows the strong chemi- 
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cal differentiation between nuclei of galaxies. Most prominent differences are observed 
between NGC253 and M82. This can be interpreted in terms of a more evolved stage 
of the nuclear starburst in M82 The chemistry of NGC4945, though claimed to be 
at an intermediate stage of evolution between NGC253 and M82, clearly resembles 
more that of NGC253 except for the lack of SiO. Similarly, the position within the 
nuclear region of IC 342, where most of the molecular observations have been made, 
show relative abundances close to those in NGC 253. This position is far from the main 
region of star formation where PDRs are claimed to dominate the chemistry. 

As far as Maffei2 and NGC 6946 are concerned, the available molecular observations 
do not yet allow us to carry out a similarly detailed comparison. 

3. A comparison of the molecular abundances of NGC 253 with selected Galactic sources 
shows a striking similarity between NGC 253 and SgrB2(0H). The latter source, lo- 
cated in the envelope of SgrB2, is commonly taken as the prototype of the molecular 
cloud complexes in the Galactic center region. This indicates that the chemistry of the 
nuclear molecular environment in NGC 253 and the Galactic center molecular clouds 
are driven by similar chemical processes. If this is the case, the chemistry and the 
heating of the nuclear molecular material in NGC 253 is dominated by large-scale low 
velocity shocks. 

4. The comparison with Galactic sources also shows a close similarity between the abun- 
dances within the nuclear starburst of M 82 and those of the Orion Bar, the prototype of 
a photon dominated region. This resemblance fully supports the idea of the chemistry 
in M 82 being dominated by PDRs. 
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A. Gaussian fitting to blended lines 

Table 1 includes a label in the third column indicating if and how the respective spectral 
feature is affected by blending. In the following the encountered different cases are described: 
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• & — > blended: Some observed features are composed of several transitions from different 
species. If the line is only partially blended, it is still possible to separate the individual 
components by fitting different Gaussian profiles. For these fits, sometimes it was 
necessary to fix the velocity and/or linewidth. 

• s — > synthetic: The blending may prevent us from making a reliable multiple Gaussian 
fitting to an observed feature. In this case, when we have observed several transitions 

of one of the species contributing to the blended feature, it is possible to use these to 
estimate the species' contribution to the blended feature by interpolating or extrapo- 
lating the linear fit to the rotation diagrams (sec Appendix B.l). This is equivalent to 
using Eq. B6 to estimate the expected integrated intensity of a line from the derived 
column density and rotation temperature. Using the estimated line intensity and the 
linewidth and radial velocity of the other observed transitions we subtract the expected 
synthetic Gaussian profile from the spectrum. The estimated integrated line intensi- 
ties using this method have been plotted in the corresponding rotation diagrams with 
open symbols (Fig. 5). The transitions from other species contributing to the blended 
feature have been fitted to the residual spectrum after the subtraction of the synthetic 
profile. 

• hf ^ hyperfine: Some species show transitions which consist of a group of unresolved 
hyperfine components. The relative intensities of these lines are fixed by their spectro- 
scopic parameters, assuming LTE and optically thin line emission. Thus, the fitting 
has been carried out with a comb of Gaussian profiles at the rest frequencies of the 
hyperfine components with the same widths and with fixed line intensities relative to 
the main component which is taken as a free parameter. The derived parameters shown 
in Table 1 correspond to the integrated intensity of all the hyperfine components but 
the frequency, line intensity and velocity refers to the main hyperfine component. 

• m — ^ multi-transition: Some features of species such as CH3OH and CH3CN consist 
of groups of overlapping transitions, involving different energy levels. In these cases, 
a single Gaussian profile has been fitted to the observed spectral feature. The total 
integrated line intensity of all the blended transitions has been used to derive rotation 
diagrams as described in Appendix B.2. 
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B. Rotation Diagrams 
B.l. Basics 

Assuming optically thin emission for the observed molecular transitions, the column 
density in the upper level can be derived as 



where A^i, the Einstein coefficient for a transition connecting the levels u and I, can be 
written as A^i = {6A7r^u^/3hc^)SiJ,li/g^, and J^{T) = {hu / k)[exp{hu / kT) - with Tbg 
denoting the temperature of the background, Tex the excitation temperature of the transition 
and Qu the degeneracy of the upper level. 

In the Ray leigh- Jeans approximation {hf <^ kT), and assuming T^^ ^ 3bg, eq. Bl can 
be reduced to 

-^-^-^ (-) 

where the integrated brightness temperature J dv is denoted as 1^. If LTE is as- 
sumed, the population distribution of all the levels can be represented by a single tempera- 
ture, the rotation temperature {Tj-ot), given by 

where is the total column density, Eu is the energy of the upper level, and Z the 
partition function calculated as 

Using Eq. B2 in Eq. B3 we obtain 

TV = fl^f^ ^e^"/^^-- (B5) 

or the equivalent logarithmic expression 



Vi 
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which constitutes the basic equation for the rotation diagrams (see Blake et al. 1987; 
Turner 1991; Goldsmith & Langer 1999, for a complete description of this method). Plotting 
Nu/Qu for each transition against the energy of its upper level in a logarithmic scale, the 
population distribution can be fitted by a straight line that provides the total column density 
divided by the partition function as well as Trot- 



B.2. Blended lines 

Whenever a molecular feature is composed of a number of unresolved blended transitions 
of a given molecular species, rotation diagrams cannot be directly used to determine the 
physical properties of the molecular emission. The spectroscopic parameters of the blended 
individual transitions within an observed feature allow us to estimate the contribution of 
each transition to the total integrated intensity, assuming LTE and optically thin emission. 
Thus, using eq. B5 we can express the relative integrated intensity between two transition 
with Einstein coefficients A'^i and A'^i, energies of the upper level E'^ and E'^, and degeneracies 

This allows us to compute the relative intensity of transition for a given T^^^. Wc will 
differentiate between the temperature, T^^t^ assume for estimating the contribution of 
each transition to the blended line, and the rotational temperature, Tj-ot, derived from the 
rotation diagram. 

If we consider the case of a molecule with hyperfine structure where all the transitions 
within a feature have the same upper level energy, then we can make the assumptions of 
E'^ ^ E'l and z/' ^ v" . Therefore, the contribution of each transition will not depend on the 
assumed T^t and will only depend on the spectroscopic parameters of each transition in the 
form 

\ _ ^Ul9u pg^) 



W" J A" n"' 

/HF ^ul9u 

When the contribution of each hyperfine transition to the total line is computed and 
plotted on the rotation diagrams, they will all lie in the same point in the rotation diagram 
(see NO and NS diagrams). 

In the general case where E'^ ^ E'^, the calculated contribution of each line for a given 
T^°t, will lie on a straight fine when included into the rotation diagram. The fine traced by 
these points will have the slope logiQe/ kT^^^.. 
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The LTE approximation assumes that all the transitions will lie in a straight line in the 
rotation diagram, thus we can modify the assumed T^°^ used to calculate the contributions 
within the blended lines to obtain the best fit to a straight line to all data in the rotation 
diagram. At this point we will obtain the situation in which T^°^ Tj-ot- 

In this paper we have several examples where this method has been applied, i.e. the 
CH3CN rotation diagram, where only two multiple-transition lines are plotted, the CH3CCH 
diagram, where three blended lines are perfectly aligned to a straight line and the CH3OH 
diagram, where a mixture of multi-transition lines and single lines are plotted altogether 
(Fig. 5). 

Wc have to take special care when dealing with the symmetric top rotors such as CH3CN 
and CH3CCH. In this case, the rotation temperature is not sufficient to describe the relative 
population of the levels given that both T^in and Tj-^t are required to evaluate the partition 
function as described by Turner (1991). Thus T^t would characterize the population of 
J levels within a K ladder, while Tkin would describe the population between X-ladders. 
Our observations do not resolve the transitions of different X-ladders. We assume both 
temperatures to be the same (Tkin = T'rot) so that the population of all levels is described 
by a single rotation temperature. Given that Trot is a lower limit to Tkin, with typically 
Trot Tkin, this assumption may result in an ovcrcstimation of the column density derived 
from rotation diagrams as the contribution to the total integrated intensity of the K > 1 
levels may have been underestimated. On the other hand, the evaluation of the partition 
function will also be underestimated. Thus, the overall effect in the total column density 
determination is not expected to be more than a factor of 2 — 3 while the derived rotation 
temperature will only be marginally affected. 

C. Notes on individual molecules 

In the following, we will discuss the fitting procedure of the identified molecular species 
in detail. 



C.l. Carbon monosulfide: CS 

Only the J — 3 — 2 transition of carbon monosulfide at 147 GHz is observed in the 
survey. It appears to be partially contaminated by the H35a recombination line. The 
estimated Gaussian profile derived from the observed H34q! (Sect. C.26) has been subtracted 
to account for its contribution. 
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Two rare isotopes, ^^CS and C^^S at 144.6 and 138.8 GHz respectively, are also observed 
in the spectral scan. Martin et al. (2005) used these observations of CS and its isotopes, as 
well as additional CS lines, to derive isotope ratios of ^^S/^^S~ 8 ± 2 and ^^S/'^^S> 9. From 
the observation of additional CS transitions at 3, 1.5 and 1.2 mm, Martin et al. (2005) derive 
a Trot ^ 10 K. 

C.2. Nitric oxide: NO 

One rotational transition of nitric oxide (150.2 and 150.5 GHZ) has been identified in our 
survey. Martin et al. (2003) confirmed this identification by observing additional transitions 
at 250 GHz. These transitions have also been included in the rotation diagram in Fig. 5. 
The topology of the NO energy levels is similar to that of NS. Transitions connecting upper 
and lower levels are labeled H+ and H^. 

It was only possible to fit the H+ hyperfine group of the transition at 150.2 GHz, given 
that the H^ group is blended with the more intense H2CO line (Sect. C.14). As shown in 
Table 1, the H~ group of transitions has been estimated by fixing width and radial velocity 
of each component to the values derived from the fitting to the H+ transition and the line 
intensities determined by their relative line strengths. 

C.3. Nitric sulfide: NS 

One rotational transition of nitric sulfide (161.2 and 161.7 GHz) has been detected in 
our survey. Martin et al. (2003) confirmed the identification by detecting other transitions 
at 115 and 207 GHz. The rotation diagram in Fig. 5 includes these additional transitions. 

This molecule presents A-doubling and hyperfine splitting. Because of the A-doubling, 
each rotational level is divided into two levels with opposite parity. Transitions connecting 
the upper levels are denoted as e and those connecting the lower levels as / (see Table 1). We 
are able to resolve the spht rotational levels due to A-doubling but the intrinsic hnewidths 
prevent us from resolving the hyperfine structure. 

C.4. Sulfur monoxide: SO 

We detect three transitions of sulfur monoxide (138.2, 159 and 172.2 GHz). The veloc- 
ity components at 180kms~^ and 275kms~^ are clearly identified. As seen in the rotation 



-58- 



diagrams (Fig. 5), the 275kms~^ component of the 84 — 23 transition shows a much higher 
intensity than expected if we assume both components to have the same rotation tempera- 
ture. This is hkely due to the contamination of this component by the Si^g ~ 7ij transition 
of NH2CN. The uncertainty in the contribution of this hne due to the noise strongly affects 
the fitting of the high velocity component of this SO transition. 

One of the isotopic substitutions of sulfur monoxide, ^^SO (135.8 GHz) is detected 
blended with a SO2 transition. This faint feature has been fitted with a double Gaussian, 
with the velocity and width of each component fixed to the values derived from the same 
transition, 4,3 — 82, of the main isotope. We calculate a ratio ^^SO/^^SO= 5.1 ± 1.2, using 
both velocity components of the A3 — 82 transition of both isotopes. This ratio is consistent 
within the errors with ^28/343= 8 ± 2 as derived by Martm et al. (2005) from CS data. 

C.5. Silicon monoxide: SiO 

Two transitions of silicon monoxide in the ground vibrational state are observed at 180.8 
and 178.7GHz (J = 8 — 2 and J = 4 — 8, respectively). It is possible to identify two velocity 
components from the J = 8 — 2 transition at 180 and 260 kms~^. Velocity and width of the 
fainter and more noisy J = 4 — 3 transition were fixed to the values obtained from the fit 
to the J = 3 — 2 line. The rotation diagram has been derived for each velocity component 
(solid and dashed lines in Fig. 5). 

C.6. Hydrogen sulfide: H2S 

Only one line of hydrogen sulfide at 168.8 GHz, with velocity components at 180 and 
275kms~^, is detected in the 2 mm survey (see also Martin et al. 2005). Physical parameters 
in Table 2 have been derived by assuming a rotation temperature of 12 K similar to that 
derived for many other species (see Table 2). 

C.7. Hydrocyanic acid: HCN 

The J = 2 — 1 transition of HCN at 177.2 GHz lies outside the covered frequency range, 
but the 2 — 1 hne of its H^^CN isotope at 172.6 GHz was detected. A fit to this transition 
allows us to estimate the abundance of HCN if we assume a rotation temperature of 12 K 
and the ^'^C/^^C ratio of 40 derived by Henkel et al. (1993). 
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C.8. Oxomethylium: HCO+ 

The only 2 mm transition of the main species of oxomethyhum, with energy levels low 
enough to be observable in the ISM of NGC 253, at a frequency of 178.3 GHz, lies a few GHz 
above the observed frequency range. Nevertheless, one of its isotopes, H^'^CO+, has been 
detected in the J — 2 — 1 transition at 173.5 GHz, near the upper end of the spectral scan. 
Two velocity components can be fitted to the observed profile. The HCO''" column density 
shown in Table 2 has been calculated from the integrated intensity of the measured H^^CO+ 
line with an excitation temperature of 12 K and an isotopic ratio of ^^C/^^C ~ 40 (Henkel 
et al. 1993). 

The J = 2 - 1 transition of HC^^O+ is also identified at 170.3 GHz. The low velocity 
component shows a higher intensity than expected from the observed H^^CG"*" profile as well 
as a slightly higher velocity. If we compare both profiles we derive a H^^GO^/HG^^O^ ratio 
of 1.7 ± 1.1 and 2.9 ± 1.2 for the low and high velocity components, respectively. The value 
derived for the high velocity component agrees within the errors with the expected value of 
~ 3.7 if we assume the isotopic ratios ^^C/^^C = 40 (Henkel et al. 1993) and ^'^O/^^O = 150 
(Harrison et al. 1999). 

C.9. Hydroiso cyanic acid: HNC 

Similar to HGN, the J = 2 — 1 transition of HNG at 181.3 GHz lies outside the frequency 
range of the survey. We identify a feature at 174.1 GHz as the J = 2 — 1 transition of HN^'^G. 
The detection is uncertain and thus tentative as the line is strongly blended with the also 
tentatively identified 3i,2 — 2i,i transition of C-G3H (Sect. G.18). 

C.IO. Ethynyl: C2H 

The J = 2 — 1 group of hyperfine transitions of G2H, observed at 174.7 GHz near the 
upper frequency cutoff of the survey, is the second brightest feature observed in the 2 mm 
window after GS. Even though the hyperfine structure is unresolved, it is still possible to 
differentiate the two velocity components at 160 and 280kms~^. 
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C.ll. Thioxoethenylidene: C2S 

The survey shows four faint hues of thioxoethenyhdene (131.6, 140.2, 142.5 and 144.2 GHz). 
The low intensity of the observed transitions makes the fitting uncertain (see Table 1). Only 
the llio — IO9 and lOn — 9io transitions at 131.6 GHz and 140.2 GHz are reliably detected, 
while the transitions at 142.5 GHz and 144.2 GHz are tentative. 



C.12. Carbon oxide sulfide: OCS 

Two transitions of carbon oxide sulfide (133.8 and 145.9 GHz) are detected in the survey. 
In the rotation diagram in Fig. 5 these observations are complemented by the 3 mm OCS 
transitions presented by Martin et al. (2005). 

The J = 13 — 12 transition of OCS at 158.1 GHz might have been detected close to a 
SO2 line, but the signal-to-noise ratio is not high enough for a reliable Gaussian fit. 

C.13. Sulfur dioxide: SO2 

Five transitions of sulfur dioxide (134.0, 135.7, 140.3, 146.6 and 151.4 GHz) have been 
identified. The upper energy levels of the transitions range from 12 to 43 K providing a good 
estimate of the rotation temperature. 

The 5i,5 - 4o,4 transition at 135.7 GHz is blended with a ^''SO line (Sect. C.4). Veloc- 
ity and width of the fine were fixed to perform a two Gaussian fit for the SO2 and ^'^SO 
transitions. 

The 22,0 — 2i,i fine at 151.4 GHz is contaminated by emission of a much fainter C-C3H2 
line which, according to its estimated intensity (Sect. C.21), would account for ~ 15% of the 
observed feature. The contribution of the C-C3H2 strongly depends on the assumed T^-ot- If 
we consider the uncertainty in the Tj-ot derived from C-C3H2 and estimate the contribution 
of this transition to the observed SO2 line, the effect on the derived SO2 parameters would 
be of ±1 K in Trot and ±5% in the column density. 

The 52,4 — 5i,5 and 7ij — 60,6 transitions are blended with a CH3OH group of lines at 
165.1GHz. We have estimated their contribution to the observed feature by assuming a 
Trot = 15 K as derived from the rotation diagram. 

For the derived excitation conditions of SO2, its 32,2 ~ 3i,3 line at 158.1 GHz is expected 
to have an intensity of ~ 5 mK, which is below the noise level at this frequency. 
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C.14. Formaldehyde: H2CO 

Two ortho (140.8 and 150.5 GHz) and one para (145.6 GHz) J = 2 — 1 transitions of 
formaldehyde are detected. In each of the profiles, the two velocity components at 180 and 
285 km s^"*^ are clearly separated. 

The 2o,2 — lo,i transition at 145.6 GHz and the 2i^i — li^o at 150.4 GHz are blended with 
an HC3N and a C-C3H2 line, respectively. Before performing the Gaussian fits to the H2CO 
features, the estimated contributions from HC3N (Sect. C.19) and C-C3H2 (Sect. C.21) were 
subtracted from the spectra. 

Rotation diagrams are plotted for each velocity component (Fig. 5) where rotation tem- 
peratures of 27 K and 34 K are derived for the 180 and 285 km s^"*^ components, respectively. 
The 285kms~^ component of the 2i^2 — li,i transition appears to have a significantly higher 
intensity than expected from the rotation diagram. The large intensity of the 2i 2 — li,i line 
of H2CO could be due to the contamination by an unidentified line. If this measurement 
is not taken into account, the resulting temperature derived for the 285kms~^ component 
would be 28 ± 2 K while the estimated column density would be lower by ~ 30%, closer to 
the parameters derived from the 180kms~^ component. 

C.15. Thioformaldehyde: H2CS 

Three transitions, all belonging to ortho-thioformaldehyde (135.3, 139.5 and 169.1 GHz), 
are tentatively detected. 

The 4i,4~3i 3 transition is blended with the H36q; recombination line. A Gaussian profile 
similar to that of the observed H34q; fine at 160 GHz (Sect. C.26) has been subtracted. We 
assume that the residual emission is due to H2CS. We also find the 4i_3— 3i_2 line blended with 
the CS J = 3 — 2 emission from the image band. With the known image sideband rejection 
at this frequency (Fig. 3), we can estimate the contribution of the CS line measured at 
146.9 GHz (Sect. C.l). The resulting feature after the subtraction of the CS line is fitted by 
a Gaussian profile. There is a part of this feature which is not properly fitted by emission 
from CS and H2CS. This may be caused by a slight difference in the observing position 
between this spectrum and that containing CS in the signal band, which would cause an 
appreciable change in the shape of the subtracted line. Contamination by a third species is 
unlikely. The 5i^5 — 4i 4 fine at 169.1 GHz is not blended, but the relatively high noise level 
at this frequency and the presence of nearby lines makes the fit uncertain. 

A rotation temperature of 11 K is determined from the rotation diagram (Fig. 5) and 
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is consistent with those derived from other molecules. With this rotation temperature, the 
4o,4 ~ 3o,3 transition of para-H2CS at 137.3 GHz is expected to have an intensity of ~ 6 mK. 
The non-detection of this line could be explained by an excitation temperature of > 20 K. 

C.16. Isocyanic acid: HNCO 

The two brightest lines of isocyanic acid in the covered frequency band are clearly 
detected (131.9 and 153.9GHz). The 280kms~^ component of the 7oj — 6o,6 transitions 
at 153.9 GHz is blended with a group of CH3CCH lines. Thus the contribution of this 
component has been calculated by assuming the 180 and 280kms~^ velocity components 
have the same excitation temperature. This assumption is supported by the good fit of 
the CH3CCH line to the residual feature (Sect. C.25). The line observed at 154.4GHz is 
tentatively identified as the 280 km s^^ velocity component of the HNCO 7i. 6 — 61,5 transition. 
The 7i J — 61,6 transition at 153.3 GHz, which should have a similar intensity, is not detected. 
Due to the uncertainty of this identification, this line has not been used in the analysis of 
the HNCO excitation conditions. Including this transition in the rotation diagram (Fig. 5) 
would result in a derived Trot ~ 29 K, slightly higher than that derived for the 180kms~^ 
velocity component, and a column density ~ 16 % lower than that in Table 2. 

C.17. Protonated carbon dioxide: 11000+ 

We tentatively identified two transition (149.6 and 171.0 GHz) of protonated carbon 
dioxide, a quasilinear molecular ion and nearly symmetric top. These are the brightest 
observable transitions in the 2 mm range. Only the high velocity component of both lines 
are detected. The large intensity observed in the 81 § — 7i.7 transitions, expected to be similar 
or lower than the 7qj — 60,6 j is mainly due to the uncertainties of the baseline substraction. 
A rotational temperature of 12 K have been assumed to estimate an average column density 
with both detected transitions. 

Further observations of other HOCO"'" transitions are needed to confirm the detection 
and constrain the excitation temperature of this molecule. 

C.18. Cyclopropenylidyne: C-C3H 

The 3i,2 — 2i_i group of transitions at 174 GHz is the only feature of C-C3H identified in 
the survey. Contamination by the feature by the J = 2 — 1 transition of HN^^C makes the 
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fitting uncertain. Thus, the velocity component were fixed to Vlsr = 180 and 280kms~^ as 
derived for most species and the width of the lower velocity component to At;i/2 = 41 kms~^ 
similar to that from the fitting to the high velocity one. The 4i 4 — 3i 3 transitions at 
172.6 GHz appear to be blended with a H^^CN fine. At 132.9 GHz, the 3i,3 - 2i,2 group of 
transitions, which is expected to be similar in intensity to those at 174 GHz, is likely present 
blended to a transition of CH2NH and an also likely transition of CH3OH, but the noise level 
at these frequencies do not allow us to separate these lines. Given the considered source size 
of 20" the 133 GHz line would be 20% fainter than that at 174 GHz, so beam dilution cannot 
account for the non-detection of these fines. Therefore, its detection, similar to the HN^^G 
line to which it is blended (Sect. G.9) is highly tentative. 



C.19. Cyanoacetylene: HC3N 

Four out of the five transitions of cyanoacetylene detected in the observed frequency 
range (136.5, 145.6, 154.7, 163.8 and 172.8 GHz) appear not to be blended with any other 
molecular line. Two velocity components have been fitted to each fine. The J = 16 — 15 
transition at 145.6 GHz is blended with an H2GO fine. Since we have observed more HG3N 
than H2CO lines to accurately determine the physical parameters, we have estimated the 
relative contribution of each of the velocity components of HC3N. For the HC3N J = 18 — 
17 line, where the velocity components arc not clearly difi^erentiated, the velocity of each 
component has been fixed. The J = 19 — 18 transition seems to lie above the T^t ~ 
20 K derived from the other observed transitions. This may indicate that the lower energy 
transitions would not be optically thin as blending with any other molecular line seems 
unfikely. Multifine observations by Mauersberger et al. (1990) (towards a position 9" north of 
our position) also noticed an excess of the higher transition lines of HG3N. Their comparison 
with model computations show that two molecular components are required to explain the 
observed line intensities, one with n(H2) ~ lO^cm"^ and Tkm ~ 60 K, and a second with 
n(H2) ~ 5 lO^cm"^ and T^in ~ 150 K. The line intensity of the J = 17—16 transition 
measured in this work is about half that measured by Mauersberger et al. (1990) which 
indicates the strong effect of a pointing offset in the total integrated intensity. However, 
as far as the physical parameters are concerned, the use of the J = 19 — 18 fine, although 
causing a change in the determined rotation temperature of ~ 50%, affects the total column 
density by less than a factor of 2. 
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C.20. Methanimine: CH2NH 

One transition of methanimine is clearly identified in the line scan. The observed profile 
of the 2i^i — li,o transition at 133 GHz is heavily affected by non linear baselines and the 
fit shown in Table 1 is uncertain. The li_o ~ lo,i transition at 167 GHz, expected to have 
an intensity twice as large as that at 133 GHz, is not rehably detected due to the lack of 
sensitivity at this frequency. Although the detection seems to be clear, further observations 
of transitions in other frequency bands are needed to confirm this detection. 

C.21. Cyclopropenylidene: C-C3H2 

Three transitions of cyclopropenylidene (155.5 and 150.8 GHz) are clearly detected. 
However, only one transition, the 32,2 — '^1,1 hne at 155.5 GHz, appears not to be blended 
with any other line. 

The 4o,4 — 3i,3 and 4i,4 — 3o,3 lines, both at 150.8 GHz, appear to be blended. A double 
Gaussian profile was fitted by assuming equal width and velocity and a line intensity ratio 
determined by their spectroscopical parameters, assuming optically thin emission. 

The intensities of the other blended transitions (145.1 and 150.4 GHz) were estimated 
for a rotation temperature of 9 K as derived from the rotation diagram (Fig. 5). Given the 
small dynamic range in energies of the unblended transitions used in the rotation diagram, 
the estimated Trot is affected by a large uncertainty (see Table 2) and, therefore, the estimates 
for the blended transitions are also quite uncertain. The 5i,4 — 5o,5 transition at 151.3 GHz 
with an estimated intensity of ~ 1 mK, well below the noise level of the spectra, has been 
calculated with the aim of estimating its contribution to the faint SO2 transition at that 
frequency (Sect. C.13). 

C.22. Cyanamide: NH2CN 

Only one line of cyanamide is clearly detected in the survey, the 81 j — 7i,6 transition at 
161.0 GHz which appears to be partially blended with image sideband emission from H2S. 
The emission of H2S in the signal band has been corrected by the image sideband rejection 
and subtracted from the spectrum before fitting the NH2CN line. The other observed lines 
are close to the noise level and therefore all fitted parameters are strongly affected by the 
baseline. The 81,8 — 7ij transition at 158.8 GHz, seems to be blended with a SO (Sect. C.4) 
but the noise level does not permit a fit. The high rotation temperature of 63 K derived 
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from the rotation diagram (Fig. 5) depends mainly on the fainter and therefore less reliably 
fitted transitions. A much lower rotation temperature would lower the column density by a 
factor of about 2. 



C.23. Methyl cyanide: CH3CN 

Two spectral features are identified as methyl cyanide transitions (147.2 and 165.6 GHz). 
Each of the J — J' transitions of this symmetric top consist of a number of overlapped K- 

componcnts {K = ... J — 1). Single Gaussian profiles were fitted to the observed lines. 
The contribution of the K = Q and 1 ladders, for the low rotation temperatures derived from 
other high dipole moment molecules, account for ~ 98% of the line profile. Therefore the 
K > 1 contribution to the total line intensity has not been taken into account. We derive 
a ^rot = 10 K from the rotation diagram in Fig. 5 which has been obtained as explained in 
Sect. B.2. 



C.24. Methanol: CH3OH 

A total of nine transitions or groups of transitions of methanol is detected in the line 
survey at 143.9, 145.1, 146.4, 157, 165 and 170 GHz. The transitions at 145.1GHz are 
blended with the C3H2 81^2 — 22,1 line whose contribution has been estimated (Sect. C.21) 
and subtracted from the observed feature. In the same way, the lines at 165.0 GHz are 
slightly blended with SO2 transitions. The emission of these SO2 hues, subtracted from the 
spectra (Sect. C.13), mainly affects the fit of the 4i^3 — 4o,4 transition. The 6_i — 5o transition 
of CH3OH expected to be detected at 132.8 GHz, is heavily affected by baseline instabilities. 
The rotation diagram derived from the observed transitions has been plotted using both the 
single- and multi-transition lines as described in Sect. B.2. The derived rotation temperature 
is 12 K. 



C.25. Methyl acetylene: CH3CCH 

Three methyl acetylene features are observed (136.7. 153.8 and 170.9GHz). As in the 
case of CH3CN each of these corresponds to a A J = 1 transition consisting of a number of 
unresolved i^'-components (with K = Q . . . J — 1). Single Gaussian profiles have been fitted 
to the observed profiles. 

The best fit to the rotation diagram as described in Sect. B.2 results in a rotation 
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temperature of ~ 62 K. Given the high rotation temperature derived only contribution of 
transitions with K = 0...4, which have upper level energies < 160 K, have been taken into 
account in the fit. For a temperature of 60 K, transitions with K = 4 contribute ~ 5% to the 
observed profiles, while higher K lines represent less than 1% of the integrated intensities. 
Neglecting the contribution of these higher energy transitions (i.e. assuming the whole 
component is dominated by the K = transitions) would result in an overestimation of the 
column densities by up to a factor of 4. 



C.26. Hydrogen recombination lines 

Four Ha recombination lines are present within the frequency range covered by our 
survey. Only one of them, H34a, is not blended and can be fitted. The H33a line lies close 
to the 2 — 1 line of C2H at the upper end of the frequency range where the rms is similar to 
the expected intensity of the fine. The H35q; and H36q; lines are blended with CS 3 — 2 and 
H2CS 4i^4 — 3i_3 fines, respectively. If the relation of the fine flux Si, oc ul appfies at these 
frequencies (Puxley et al. 1997), the lines should have flux densities 20 and 30% below that 
measured for the H34ci; line. We have considered them to have a similar line proflle as that 
of the measured H34q; line given the uncertainty in the fitted intensity. 

To convert the measured H34q; line temperature into flux densities we use the conversion 
factor for the 30 m telescope S/Tmb — 4.95 Jy/K. Following Puxley et al. (1997), a source 
to beam coupling factor {0l + 9l)/dl has been applied for a source size ds = 7". We calculate 
an integrated line fiux of (1.08 ± 0.08) 10~^^ Wm~^. As seen in Fig. 9, this value closely 
follows the best fit to the relation between the integrated line fiux versus frequency derived 
by Puxley et al. (1997). In this figure the fiux of the H34a line is plotted together with the 
recombination lines observed in millimeter and centimeter wavelengths (Puxley et al. 1997, 
and references therein) . 



io' 

Frequency (GHz) 



10^ 



Fig. 9. — Recombination line fluxes of measured Ha lines in the range from 1 to 160 GHz as 
a function of frequency. A best-flt to the line fluxes in the form 5'l oc is shown as fltted 
by Puxley et al. (1997). 
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